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'LOC'-vf STRUCTURES IN AMORPHOUS Ge Se, SYSTEM 1 

X 1 — X 


In this thesis, the results of the investigation 
undertaken to record and analyse the X-ray absorption 
(K edge EXAFS), fluorescence and diffraction data of 
amorphous materials (Ge x Se^_ x s y si::am with x = 0.1, 0.3 
and t„0) are reported., These materials belong to the 
erno ’nhous semiconductor class of non— crystalline matter. 

Extended X-ray absorption fine structure (EX\FS) 
refers to the oscillations of X-ray absorption coefficien 
on fh.' high energy side of the absorption edge. Interest 
in i ' •■’S has been revived in past several years for its 
applications to the structure determination of amorphous 
materials. M though th^ revival of interest in EXAFS is 
due to the improvements in theory and experimentation, 
in particular the synchrotron related developments, we 
have in our laboratory recorded the EXiFS data with a 
conventional X-ray source. By operating the tube at 



higher ratings the requirement of high photon flux was 
mot. Adequate resolution (~ 4 eV) in the vicinity of 
absorption edge was obtained by using an exit slit of 
0.05 mm width (which allows minimum increment of 0.01° 
in 2 Q ) . 

The analysis of F.XAFS data involved number of 
steps with extensive use of computer. The programmes 
nw ded for the analysis have been developed. The experi- 
mental data on crystalline and amorphous germanium served 
to standardize the analysis procedure. The study of amor- 
phous <Je,/3e^ system has been limited to compounds with 
low concentration of germanium. EXAFS is particularly 
useful in such a system as we can probe the environment 
of hot! kind of atoms. 

Chapter T contains some introductory remarks on 
th.' local structure in non-crystalline materials and 
met' ods for its determination. The advantages and dis- 
advantages of the relatively new technique - EXAFS over 
tlv= traditional methods; X-ray diffraction (R. D.F. measure- 
ment) in particular have been briefly outlined. The possi- 
bility ot fscertaining the kind of backscattering atom, 
particularly in system like GeSe using EXAFS has been 
examined. 


Chapter II gives an account of the experimental 
details of all the three types of measurements namely 



the X-ray absorption (SXAFS), fluorescence and diffrac- 
tion. Dispersion and resolution of the spectrometer has 
also been discussed. 

In chapter III the theory underlying the SXAFS 
phenomenon and data analysis procedure has been presented 
bri ?f ly. Chapter IV presents the results obtained by 
analysing the experimental data of germanium. K edge 
SXAFS of both crystalline and amorphous germanium have 
been analysed. The effect of width of the filtering 
window, used to separate out the first shell EXAFS, on 
the phase and amplitude functions and subsequently on 
the determiniation of structural parameters has been 
discussed. The phase function was found to be least 
affected by the width of the window function as indi- 
cated by the same value of near neighbour di-s-tanrcxss 
obtained with different window widths. On the other 
hand the amplitude function was found to be sensitive 
to the changes in window width. The concepts of chemical 
transferability of phase and amplitudes have been used 
to calculate the near neighbour distance and the coor- 
dination number for amorphous germanium. The determination 
of near neighbour distance by comparison of experimental 
and theoretical phase shifts (Teo, B.K. and Lee, P.A. 
(1979), J. Am. Chem. Soc. 101 , 2815) required that the 
phase functions be obtained using critical value of 
E q - E c at which the phaso shift is linear and the posi- 
tion of the peak in the magnitude of Fourier transform 



is independent of the weighting factor (Stearns, M.B. 
(1982), Phys. Rev. S2_5, 2382). The analysis of X-ray 
diffraction data is briefly presented in the last section. 

In chapter V, the analysis of experimental data 
of amorphous Ge Se. system (x =- b.l, 0.3) is presented. 

1< edge SX1FS of both kind of atoms (Ge and Se) has been 
analysed. vThile analysing SXAFS data, the phase functions 
have been extracted (at S = E ) following the procedure 
in chapter IV. It has been assumed that the first coor- 
dination shell around central atom can have both kind of 
atoms on account of the proximity of atomic numbers of 
Ge and Se. For determining the near neighbour distances 
two different cases are considered depending upon the kind 
of a tom as backscattering atom, the result obtained using 
the experimental phase function could be interpreted in 
terms of the probability of that atom being the nearest 
neighbour to the central atom. The calculated probabili- 
ties have been used in the proposed model which is based 
on the probabilities of all the different possible confi- 
gurations of the first coordination shell around the 
central atom. The experimental first shell EXAFS then 
coul-5 be compared with the one obtained from the proposed 
model to determine o set of structural parameters. DRDF 
analysis has been briefly presented in the last section. 

Chapter VI gives the summary of the present work 
and also the scope for future work. 



CHAPTER I 


INTRODUCTION 
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1.1 Non Crystalline Materials 

Solid materials can be broadly classified into two 
categories depending on their structure. The crystalline 
materials have a three dimensional periodic arrangement 
of their constituent atoms or in other words long range 
order (lro) exists. Structurally amorphous, non-crystal- 
line, disordered, glassy and vitreous are synonyms, but 
the last two terms are generally reserved for non— crystal- 
line solids prepared by quenching from liquid state. The 
amorphous materials lack the regularity of arrangement of 
atoms as in crystals. Diffraction experiments (X-rays, 
Neutron, Electron) show that the order within the amorphous 
phase does show some sort of similarity with that of the 
crystal lattice, losing however, the long range order. 
Amorphous solids distinctly differ from crystalline one 
due to lack of long range order and hence possess no 
readily definable translational symmetry. Figure 1.1 
illustrates two different phases : crystalline and amorphous 
of the same chemical composition A 2 O 2 . There is no sharp 
dividing line between the crystalline and the so called 
amorphous materials. As the crystalline size diminishes 
the widths of the lines composing the powder X-ray diff- 
raction pattern increase until the fainter lines are in- 
discemable and the stronger assume diffuse characters 
which is typical of non-crystalline substances. The 
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crystallite dimensions calculated from such broad bands 
are often of the order of one or two unit cell dimensions, 
a result very difficult to interpret on the basis of 
repetitive structure pattern. 

Amorphous semiconductors are well known examples 
of highly disordered materials. There has been in recent 
years a surge of interest in them, particularly after the 
publication of Ovshinsky (1968) describing the various 
types of switching phenomena which can be observed with 
a large number of amorphous solids. The technological 
importance soon became evident. Applications of amorphous 
semiconductors either realized or anticipated include a 
very wide spectrum (Feinleib et al. 1971? Hilton et al., 

1 966a; Ovshinsky and Fritzsche, 1973; Neale and Aseltino, 
1973? Dessaner and Clarks, 1965) such as optical recording 
devices, optical lenses, windows and filters, continuous 
dynode electron multipliers ( channel tron ) , optical mass 
memory, phase contrast holograms, high energy particle 
detectors, ultrasonic delay lines and microfiche trans- 
parencies, in addition to the exciting switching and 
memory devices and the well established applications in 
xerography. This is indeed a very impressive list, 
specially in view of the fact that the understanding of 
the amorphous semiconductors today is still far from 
satisfactory. Further the advantage of producing amorphous 
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semiconducting materials in the form of thin films makes 
them adaptable to integration with other solid state 
technologies. 

The field of amorphous semiconductors as a whole 
is of fundamental interest. Ordered materials (crystal- 
line) are known to have universal features in their 
electronic structures. It can be expected that amorphous 
semiconductors may also possess such universal features. 
Considerable efforts are being made to obtain the rigorous 
theoretical interpretations of the experimental observa- 
tions. An understanding of the properties of amorphous 
semiconductors requires not only familiarity with solid 
state physics but rather a detailed knowledge of the 
results obtained from the combined fields of physics, 
chemistry, metallurgy and electrical engineering. Consi- 
derable exploratory work on amorphous semiconductors was 
carried out during the early fifties. In recent years, 
intensive research work in several countries has been 
directed to improve the understanding of structural, 
electrical, optical, mechanical, magnetic and other pro- 
perties of these materials in order to attain the same level 
of understanding as in case of crystalline semiconductors. 

In the present work we are concerned with the structural 
aspects. 
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1.2 Determination of Local Structure 

A complete description of the properties of any 
compound requires a full knowledge of its structure. 

There are basically two main approaches for determination 
of local structure in non crystalline materials, which 
have have been used in the past, namely the diffraction 
methods and the spectroscopic methods. In the first 
group,, the interaction of radiation with matter leads to 
information which can be analysed as such in the Fourier 
space or transformed to real space, either mathematically 
as in case of X-ray or neutron diffraction or physically 
as in case of electron microscopy. 

The electron or X-ray diffraction pattern of amor- 
phous materials consist of a few broad halos rather than 
sharp Bragg reflections. A small crystalline size, strains 
and imperfections broaden the sharp crystallite reflections 
and as the crystal size gets smaller a point is reached 
where the crystallite reflections overlap and diffraction 
pattern appears like that of amorphous materials (Germer * 
and White, 1941; Piggott, 1966). The probability of find- 
ing an atom on a sphere of radius r or the Radial Distribu- 
tion Function (RDF) 4 tc r Sj’(r), ? (r) being the atomic 
density, is used to give a one dimensional description 
of the atomic distribution. When RDF is normalised. 
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it contains a number of peaks at some values of r and area 
under each peak corresponds to the number of atoms contained 
in the coordination shell. To illustrate/ figure 1.2 shows 
the RDF of Ge ^2 Se 58 as obtained from X-ray diffraction 
data by Bienenstock et al. (1976). It can be noted that 
the first and second neighbour peaks are well defined in 
this case, but that* little information is obtainable about 
more distant coordinations. 

As we know, the diffraction methods give us only a 
one dimensional information - the result of RDF and intri- 
cate separation methods using different radiations, it is 
only by means of modelling that one can obtain a three 
dimensional information. 

Zacharaisen in 1932 (Zacharaisen, 1932) in a paper 
which led to the foundations of our understanding of glassy 
structure insisted that "It must be frankly admitted that 
we know practically nothing about the atomic arrangement 
in glasses". Knowledge about the non-crystalline materials 
remained almost static in the following thirty years. 
Zacharaisen' s random network - the father of subsequent 
models of amorphous solids - remained a mathematically 
unhandable concept. Diffraction experiments of Warren 
(Warren et al., 1936, Biscoe and Warren 1938) on vitreous 
silica and boron oxide were consistent with random network 
model and similar concepts were used by Bernal in probing 
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the structures of water and simple liquids (Bernal and 
Fowler, 1933, Bernal, 1937, 1959), but the concept remained 
underdeveloped and largely unappreciated. 

There has been a rapid growth of interest in amor- 
phous solids in the past decade ranging from covalent 
glasses such as vitreous silica through the related amor- 
phous tetrahedral systems of Si, Ge and H^O. Depending 
upon the method of preparation and the nature of molecules 
concerned, one might expect variety of amorphous structures. 
Despite this variety, all amorphous substances can, follow- 
ing Zacharaisen, be treated within the same conceptual 
framework: an irregular topology of linked subunit molecule. 
Although distinctivley microcrystalline materials can be 
prepared (Wagner et al., 1968), diffraction experiments 
on the amorphous solids cannot convincingly be explained 
in lattice terms (Cargill, 1970; Mozzi and Warren, 1969; 
Moss and Graczyk, 1969). As a first approximation one 
tries to develop an idealised model which describes the 
essential structure of group of systems. Modelling two 
extreme types of amorphous structures based upon an 
idealised tetrahedral geometry and upon spherically sym- 
metrical atoms cover many of the systems of most practical 
interest. In each case, real physical models were built 
in laboratory. Later, digital computers were used to 
build and refine the models. In the present investigation 
models based on tetrahedral geometry are of interest. 
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The idealised tetrahedral network is topologically 
similar to other non-crystalline aggregates. The ground 
work for realising Zacharaisen' s random network concept 
was led in 1960's, when a series of tetrahedral models 
were built in an attempt to understand the structures of 
vitreous silica and water (Ordway, 1964; Evans and King, 
1966). The vitreous structural models may be classified 
as microcrystallite, amorphous clusters, continuous random 
network (CRN), the chain model (Hilton et al„, 19 66b) 
suggested for GeTe and chain ' crossing model (Fawcett et al.., 
1972) suggested for Ge Se 1 system (xk 0.3) and so on. 

In its simplest form, the microcrystallite structure of 
the tetrahedral ly coordinated as well as of other types 
of solids is considered to be an assembly of crystallites 
with their orientations and shapes distributed as they 
might be in recrystallised solids. The amorphous cluster 
models were developed specially for application to tetra- 
hedrally coordinated systems. In an "ideal" random network 
structure, every atom is bonded to the number of nearest 
neighbours required to satisfy the chemical valence so 
that there are no internal "dangling" bonds. The structure 
can be enlarged indefinitely without breaking bonds or 
increasing deviation of the bond angles from their ideal 
values. For example this model seems to describe adequately 
the RDF of amorphous germanium (Grigorovici and Manila, 

1969; Polk, 1971; Turnbull and Polk, 1972). 
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The model describing the RDF gives us the nearest 
neighbour separations,, the coordination numbers and the 
bond angles, all of which are vital parameters for any 
basic theory that hopes to explain in. detail the electronic 
properties of non-crystalline materials. The problem 
often encountered with this approach is that quite 
different reasonable structural models yield almost 
exactly the same first neighbour peak area so that they 
can not be distinguished. This sort of ambiguity pervades 
much of the chalcogenide glass RDF work. It is quite 
extreme in Ge-chalcogenide because the structures of 
crystalline and amorphous Ge-monochalcogenides are strik- 
ingly different ( Bienenstock, 1973a) so that the knowledge 
of former gives little guidance concerning the nature of 
the later. Ge-Se system poses the worst problem because 
the two atomic numbers are almost identical. From 
figure 1.2 it is apparent that the first peak may have 
contributions from Ge-Se, Ge-Ge and Se-Se first neighbour 
pairs. Herein lies the essence of the ambiguity in its 
interpretation. 

All the diffraction techniques directly detect 
the lro and in particular the unit of periodicity. For 
the case of non-crystalline materials, the diffraction 
techniques, as has been pointed out earlier, have only 
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limited applicability. On the other hand the other group 
of methods namely the. spectroscopic methods could be used 
as filters to select from the vast number of models those 
which have a higher probability of being real and reject^ 
those which are not compatible with, experimental results. 
The spectroscopic methods include Raman spectroscopy, NMR, 
IR and the relatively new methods such as positron annihi- 
lation and specially EXAES which provides new ways of 
probing the environment of atoms in the structure. It is 
this technique which has been used in the present work. 

1.3 Some Special Advantages of X-ray Absorption 
Spectroscopy 

Although there exist several theoretical and experi- 
mental difficulties in satisfactorily interpreting X-ray 
absorption spectra, this technique has certain special 
advantages (Cauchois, 1948) which are stated below: 

a) It is very specific regarding the element or 
the ion under investigation. 

b) All types of substances, solids (crystalline, 
amorphous, metallic and non-metallic, elemental and 
compounds), liquids and gases in different physico-chemical 
states can be studied by this method. 

c) Wide energy range is covered in the same spectra. 
With a high resolution curved crystal spectrograph, the 



13 


different features of a discontinuity such as position/ 
shape, width and fine structure give information about 
different single particle states as well as collective 
states simultaneously (Mande and Nigavekar, 1967). 

d) Information regarding the empty states in the 
conduction band is obtained in the absorption spectra 
making use of relatively narrow inner states as a scanning 
probe. 

e) It is a non-destructive technique and compari- 
tively a small amount of the material is required for the 
inves tigation. 

f) X-ray absorption spectra are sensitive to 
chemical combination and can therefore be used in the 
study of chemical bonding. 

g) Amongst the different techniques used in the 
study of the solid state, it is relatively simple and 
inexpensive. 

h) Oscillations in the absorption coefficient on 

the high energy side of the absorption discontinuity (EXAFS) 
provides information about local structure around the 
absorbing atom. 

1.4 Extended X-ray Absorption Fine Structure (EXAFS) 

The term extended X-ray absorption fine structure 
refers to the oscillations of the X-ray absorption coeffi- 
cient on the high energy side of the absoiption edge 
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extending to about 1000 eV above the edge. There has 
been, in recent years, a revival of interest in EXAFS, 

Several sophisticated theories (Stern, 1974; Lee and 
Pendry, 1975; Lee and Beni, 1977) have been putforth to 
explain the occurrence of EXAFS. Experimental methods 
for recording the X-ray absorption spectra have also 
improved in the last decade. 

In late 1930s, X-ray absorption spectra were taken 
in transmission with low power conventional tubes, a dis- 
persive Bragg spectrometer and photographic film as 
detector (Kossel, 1920; Kronig, 1931, 1932). Lytle (1966) 
improved the measurements using modern counting and auto- 
mation techniques. The replacement of conventional X~ray 
tubes with synchrotron radiation signalled the beginning 
of new era in measurement techniques (Kincaid and Eisenberger, 
1975). A transmission technique using a channel cut mono- 
chromator and an ion chamber was developed, followed by 
development of fluorescence detection technique based on 
solid state detector (Jaklevic et al., 1977). A new 
graphite array system was developed (Hastings et al., 1979) 
which extended the limits of EXAFS experiments. The 
synchrotron related developments have been responsible 
for the renewed interest in conventional rotating anode 
X-ray tubes coupled with curved crystal monochromator 
to perform laboratory based transmission experiments 
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(Knapp, 1978). Even the conventional sources coupled 
with flat crystal monochromator yield data of fairly good 
quality (Lytle et al„, 1975) and in the present work we 
have used a similar set up for recording the EXAFS data. 

The requirement of high photon flux is met by operating 
the tube at higher ratings. 

It is believed that study of EXAFS can give valuable 
and very precise information about the interatomic distances 
in all kinds of materials. The RDF studies measure distances 
between all possible pairs in a multicomponent system 
whereas EXAFS can provide information about the immediate 
environment of an individual species in the system. EXAFS 
being basically an sro phenomenon, it is most useful for 
studying the disordered materials. Recently EXAFS is being 
used in disordered systems (Crozier and Seary, 1981; 

Crozier et al., 1977), Catalysts (Lytle et al., 1980), 
Biological systems (Shulman et al,, 1975, Sayers et al., 

1976) and so on. The distances could be determined to an 
accuracy of ± 0.01 A and the coordination numbers within 
2 0% accuracy. The accuracy is mainly deterimined by the 
quality of experimental data and the data analysis proce- 
dure used to extract the relevant parameters. The main 
problem in the analysis is to obtain the smooth free 
atom absorption, J^ Q ( K), required for normalising the 
EXAFS modulations to per atom scale (Criteria for 
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automatic background removal and normalisation has been 
recently suggested by Cook and Sayers (1981 )). Any error 
in normalisation directly affects the amplitude of the 
normalised EXAES and subsequently the determination of 
coordination number. Other factors which could affect 
the EXAFS amplitudes are thickness of the absorber (Stern 
and Kim, 1981), many body effects (Stern et al., 1980) and 
the temperature. The temperature dependent Debye-Waller 
type factor is a measure of both the thermal as well as 
static disorder in the material. This factor must be known 
for determination of coordination number. In many cases the 
effective Debye-Waller approximation is inadequate in dis- 
ordered systems (Eisenberger and Brown, 1979; Crozier and 
Seary, 1980) and one has to resort to modelling for obtain- 
ing meaningful results. 

1.5 Scope of the Present Work 

Considerable attention has been provoked to the 
study of amorphous chalcogenides, specially amorphous 
Go Se. system (ifc. Tronc et al. , 1977 and Yoshimura et al. , 
1978) because of their possibilities for numerous applica- 
tions in areas such as optical storage, imaging, replica- 
tion of patterns and phase holography (Goldschmidt et al., 
1977; Zembutsu et al., 1975; Terao et al., 1972), An 
understanding of the properties of amorphous chalcogenides 
requires a full knowledge of their structure. However, the 
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local structures of these chalcogenides have not been well 
established. Attempts have been made for determining their 
local structure (Bal and Chamberlain, 1978; Bletskan 
Gerasimenko, 1978; Kawamura et al., 1982). The Ge chalco- 
genides, in particular the Ge Se- system poses difficulty 
as the bond lengths Ge-Ge, Ge-Se and Se-Se are not very 
different and hence could not be resolved (Fawcett et al,, 
1972), One could imagine, quite readily, that such a glass 
is phase separated into GeSe 2 and Ge Se. The GeSe 2 would 
have four fold coordination of Ge by Se atoms and two fold 
coordination of the Se atoms by Ge atoms, as in amorphous 
Si0 2 . The GeSe might have each atom coordinated by three 
of the other species, as discussed by Bienenstock (1973b), 
Alternatively, it might phase separate into GeSe 2 plus Ge, 
so that all Ge atoms are four fold coordinated, instead of 
some being three, and all the Se are twofold coordinated 
by Ge. Another model considers the system as homogeneous 
molecular Ge 2 Se 3 . Here, each Ge is coordinated by one Ge 
and three Se's while each Se is coordinated by one Ge and 
one Se. For Se rich alloys (x< 0.3) we can expect amorphous 
Se structure to be maintained. The structural element of 

Se lattice is an isosceles triangle (Richter, 1972). It 

o 

is determined by shortest atomic distances r^ =2.32 A 
o 

and r ^ = 3.69 A within the screw chain of lattice. The 

o 

shortest chain distance is r = 3.46 A. In one of the 
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models suggested for Se rich alloys, this Se chain struc- 
ture is maintained and Ge atoms act as chain crossing 
points (CCM model). In the other model (RCM) Ge atoms 
are allowed to bond to other Ge atoms as readily as the 
two fold coordinated Se atoms. Both of these models 
adequately describe the RDF (Fawcett et al., 1972). The 
diffraction methods thus met with several limitations. 

Now that EXAFS has emerged out as a powerful tool for 
the determination of local structures in amorphous systems 
(Gurman, 1982), it is therefore worthwhile to study the 
Ge Se, system using EXAFS. 

X X 

In this thesis, we report our investigations on 

a) EXAFS associated with K X-ray absorption discontinuity 
in crystalline and amorphous germanium. 

b) EXAFS associated with the K X-ray absorption disconti- 
nuity of germanium in Ge Se. system (x = 0.1, 0.3). 

c) EXAFS associated with K X-ray absorption discontinuity 
of selenium in amorphous Ge Se. system (x = 0.1, 0.3). 

X X *“X 

d) X-ray diffraction studies of amorphous Ge and 
Ge x Se l~x (x = 0.1, 0.3). 

e) X-ray fluorescence studies of these materials. 

Many theoretical difficulties are involved in 
analysing EXAFS data. Inspite of these difficulties the 
author has made a successful attempt to develop computer 
programs for the analysis of EXAFS data. 



CHAPTER II 


EXPERIMENTAL 
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2.1 Specimen Preparation 

The amorphous chemical compounds Ge Q g and 

Ge^ 3 S eg q were prepared as per the following procedure. 
Appropriate amounts of constituent elements were mixed and 
sealed in a quartz tube under high vacuum. The elements 
used were of five nine purity. The bulk samples were 
obtained by melting,, annealing and then quenching in air. 
Amorphous germanium was prepared/ in the form of thin film, 
by evaporating five nine pure germanium in 10 _ Torr vacuum 
on aluminium substrate. The substrate temperature was 300°K. 

The X-ray diffraction patterns of all the samples were 
recorded. The patterns showing a broad peak vanishing 
rapidly at high diffraction angle are typical of amorphous 
compounds. Figure 2.1 shows such a pattern obtained for 
Ge Q 2 Se o 9 w hich peaks about 15.5° (29). Application of 
Debye formula (Klug and Alexander, 1954) for this pattern 
reveals that the average crystalline size in this sample 
is around 11 The compositions of the constituent elements 
were confirmed by X-ray fluorescence analysis (sec. 2.5). 

2.2 Experimental Procedure - EXAFS 

The basic experiment involves the measurement of 
absorption coefficient , , as the X-ray photon energy is 

varied across, and beyond the K-absorption edge. The 
absorption coefficient is related to the logarithmic 
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ratio of the incident and transmitted intensities at a 
particular energy (or wavelength). Thus the experiment 
consists of basically the measurement of X-ray intensi- 
ties as a function of photon energy. The following 
sections briefly describe the experimental procedure for 
EXAFS measurements. 

2,2.1 Dispersion 

Since X-ray absorption experiment involves detection 
of X-rays having different wavelengths# single crystals 
are invariably employed for separating out a particular 
wavelength. Crystal dispersion arrangements may be broadly 
classified as flat-crystal and curved crystal. Actually, 
curved crystals do not focus X-rays in the true optical 
sense, but diffract them in such a way that they converge 
to a line or a point. In our experimental set-up a flat- 
crystal arrangement has been used. A single crystal of 
LiP cut at its (200) plane (2d = 4.028 A) served as an 
analyser. X-rays incident at a glancing angle ©, only a 
particular wavelength given by Bragg's formula 

2d sin© = nX (2.1) 

are detected by the detector. The harmonic contribution 
would also be there, but being negligibly small compared 
to the first order, may be neglected. 
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2.2.2. Geometry of the experiment 

Figure 2.2 shows the geometry of X-ray diffraction 
from the analysing crystal. Crystal is set in position A. 

For every wavelength a cone of divergent rays from the 
absorber intersects the crystal surface at range of angles. 
One of the rays is likely to intersect the crystal at the 
Bragg angle, as shown for the lower ray in the figure. 

The other ray is diffracted at the same point if the 
crystal is rocked. Each wavelength present undergoes the 
same phenomena at its own Bragg angle, as is diffracted 
to a different point on the spectrometer circle, resulting 
in the dispersion of the spectrum on the basis of wavelength. 
This flat crystal arrangement is called Bragg-Brentano 
parafocussing (Brentano and Landasy, 1954). 

2.2.3 Sp ec t ral ran ge 

In principle, the largest wavelength which a crystal 

can diffract is egual to twice the interplaner spacing (2d) 

of the crystal planes parallel to its surface. For LiF (200) 
, o 

this limit is ~4.02 A. The K-absorption wavelength corres- 
ponding to germanium (1.11652 A) and selenium (0.97977 %) 
lie well within this range. Knowing the absorption edges, 
the position 2© on the goniometer circle is immediately 
ascertained by Bragg relation (2.1). The corresponding 
energy (E ) of the photon is given by 

= 12398. 52/ A (A) 


E (eV) 


(3.2 ) 
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The spectral range of interest for the K-absorption 
edge is about 1000 eV on high energy side and about 200 eV 
below on the low energy side of the edge. The corresponding 
ivavelengths and Bragg angles can be calculated using rela- 
tions (2.1) and (2.2). The range on 20 scale for the K- 
absorption edge of germanium comes out to be 29.4° to 33.8°. 
The corresponding range for selenium K-absorption edge is 
28.1° to 29.3°. 

2.2.4 Preliminary adjustments 

In the present investigation, a Seifert's X-ray 
spectrometer with auto-step scanning mechanism has been 
employed. At the centre of the goniometer of radius 17.5 orris 
is a stage which could be rotated, the axis of rotation 

V 

being also the goniometer axis. The LiF crystal, after 
polishing, is mounted on this stage. X-rays from the 
scaled X-ray tube pass through entrance slit and are incident 
on the analysing crystal at Bragg angle 0. The width of the 
entrance slit can be adjusted to a desired size. A lead 
absorber was kept before the detector to avoid damage of 
detector when goniometer is in 0° position. The crystal 
holder at the centre was rotated slightly in such a way 
that maximum intensity is received at the detector. In 
this position half the beam is obstructed by the crystal 
and remaining passed parallel to its surface to the detector. 



After this 0° adjustment, the spectrometer is said to be 
aligned. 

The target of the sealed X-ray tube usually contain 
some impurities which produce emission lines at their 
respective Bragg angles. These impurity lines were used 
to redetermine the interplaner spacing (2d value) of the 
crystal. With the present adjustments the 2d value for 
LiP crystal using FeK°< line turns out to be equal to 
4.02 62 A. 

The X-ray beam, diffracted from the crystal, passes 
through another slit called the exit slit before it reaches 
the detector. The size of this exit slit determines smallest 
increment in 29 and consequently the resolution. Resolution 
is different from the angular dispersion. The resolution 
or resolving power is defined as *-/dX or E/ £>E. 

Prom (2.1) 

2dcos9 4^0 = n hA (2,3) 

Using (2.2) and (2.3) it can easily be shown that 

U>E = ..EcotQ 0 9 (2.4) 

Thus for higher resolution, intensities should be 
measured at smallest possible increment £*9 and according 
to Rayleigh's criterion©! resolution, the width of the 
exit-slit should have maximum 2 4© divergence. The 



27 


requirement of high photon flux in EXAFS experiment 
constraints one from using an exit slit of very small 
divergence. An exit slit of .05 mm width, which corres- 
ponds approximately to 0.02° divergence, was put before 
the detector. This allows a minimum increment in 2© 
equal to 0,01°. In case of K-edge of germanium, for 

instance, 2© = 32.1° and E = 11104 eV. Equation (2.4) 

o 

gives ^E -w3.35 eV. 

2.2.5 Preparation of absorber 

Because of the high absorption by all materials 
to X-rays having energies close to their absorption edges, 
the chief problem in experimental X-ray absorption spec- 
troscopy is the preparation of homogeneous samples having 
requisite uniform thickness, typically of the order of 
few microns. Ideally nothing but the substance itself 
should bo placed in the path of X-ray beam. A convenient 
method is to prepare the absorber in the form of uniform 
thin film. In case of amorphous compounds Ge^ ^SeQ g 
and Ge^ 3 S Gq q absorbers were prepared by spreading a 
fine powder of the sample on an adhesive cellotape. The 
homogenity of the absorber was tested by recording the 
transmitted intensity through different parts of the 
absorber (rotating the absorber). The fact that the 
thickness of the absorber affects the fine structure has 
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been appreciated for quite some time. Since counter has 
been employed for measuring intensities, additional cri~ 
terion is enforced by counting statistics. Assuming that 
the statistical precision of all the eixperimental points 
is the same, it was suggested (Parratt et al. 1957) that 
the thickness should be such that I/I 0 '~ 0.1 where I and I Q 
are transmitted and incident intensities respectively. 
Thickness effect on EXAFS amplitudes has been recently 
discussed by Stern and Kim (1981). Significant distor- 
tion in EXAFS amplitudes occurs when &^-i Q x ^1.5, where 
is the K-edge step in the absorption coefficient. 

In practice, however, the thickness of the absorber can 
not be controlled and several thicknesses of samples were 
tried before the final measurement was made, 

2 , 2 . 6 Recording absorption spectra 

After the preliminary adjustments discussed in 
section 2.3.4, the absorber was placed in the path of 
X-ray beam. It is generally believed that the absorber 
can be placed anywhere in the path of the X-ray beam, 
actually it is better to place the absorber ahead of the 
crystal or near the detector so that any fluorescent 
radiation emitted by the absorber does not reach the 
detector. Both the positions were tried. Since the 
transmitted intensity from the absorber depends on the 
intensity of the beam incident on it, the placement of 
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absorber ahead of the crystal reduces the transmitted 
intensity. Because of the requirement of high photon 
flux, the absorber was finally placed near the source. 

The spectral region to be scanned for a particular 
K-absorption edge was calculated (section 2.3.3). The 
detector was moved to the initial 2 Q position. The 
detector is either a proportional or a scintillation 
counter. Since the fluorescence is generally completed 
in 10" seconds, the recovery time of scintillation counter 
is excellent. Counter signals are proportional to the 
energy of X-ray photons in both counters. Scintillation 
counter was preferred because the maximum useful count 
rate, without substantial coincidence loss, for scintil- 
lation counter (10 counts /sec.) is higher than that for 
a proportional counter (10^ counts/sec.). 

The scaled X-ray tube with Cu-target was operated 

at 50 KV, 30 mA to get the maximum possible intensity. 

The intensity could be measured in either preset count 

mode or preset time mode. To have equal statistical 

precision for all experimental points, the intensities 

5 

were measured in preset count mode (10 counts). The 
goniometer was operated in auto-step scanning mode. The 
auto-reset and auto-print switches were in 'ON 1 position 
and step-scanning was activated. The selector switch was 
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kept on 'external-pulse' position. The step in 20 could 
be adjusted within a range of 0.01° to 0.99°. Since the 
exit-slit of .05 mm was used at the detector, the lowest 
step of .01° in 20 was adjusted. The print-pulse activates 
the scanning motor, and the detector is advanced to next 
20 position with a specified step. The detector is locked 
in this position till the next print command is executed, 
giving the time for preset number of counts. The time T 
with absorber in the path of the beam was printed on a 
machine role. The scan was repeated without the absorber 
to give T q at each 20 position. The intensity I Q may 
change with time, however the change was found to be 
within the statistical precision. 

Knowledge of the start position and step size in 
2 0 allows calculation of X-ray ’wavelength for each data 
pair (T, T ) and the absorption can be calculated. The 
total experimental data consisted of 350-450 data pairs. 

An elementary computer program was developed which 

incorporated i) calculation of the absorption Jo x according 

3 

to equation $ . 1 , wavelength and energy at each point, 

3 4 

ii) fitting a Victoreen curve of the form A A B X 
to absorption data below the edge and iii) extrapolation 
and subtraction of Victoreen curve from data above the 
edge. The resulting data gives the contribution to 
absorption exclusively due to K-shell. 
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2.3 Experimental Details - X-ray Diffraction 

For the calculation of Radial Distribution Function 
(RDF), the experimental data consists of X-ray diffraction 
intensity from the sample normalized to electron units. 
Normally the powdered sample is prepared in a slid cavity 
of perspex which contributes to the measured X-ray diffrac- 
tion intensity. For this reason experiment is performed 
in two steps involving measurement of intensity from 
substrate and sample mounted on substrate (slid cavity) 
separately. Following sections describe the experimental 
procedure briefly. 

2.3.1 Geometry of the experiment 

Geometry for recording the X-ray diffraction pattern 
is very similar to the one used in EXAFS experiment. Mono- 
chromatic X-ray beam, from X-ray tube with cu-target is 
obtained using a Ni filter. After passing through an 
entrance slit, beam is incident on the powder sample at 
glancing angle 9 and the diffracted beam is received at 
the detector at an angle 29 . Figure 2.3 illustrates 
schematically the geometry of the counter diffractometer. 

The sample 0 is the flat-faced briquet of powder. 
Monochromatic radiation diverges from the entrance slit E, 
which is just a line focus on the target. The diffracted 
radiation is selected by a narrow receiving slit R. Distances 
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EO and OR are equal and since the sample face is maintained 
symmetrical with respect to the primary and diffracted beams, 
there is an approximate focussing condition. The radiation 
passed by the receiving slit (or the exit slit) is detected 
by the scintillation counter D, For recording the diffrac- 
tion pattern the counter turns at an angular velocity <jO and 
the sample turns at rate o>/2 to maintain the focussing 
conditions. 

2.3.2 Correction factors 

The X-ray diffraction intensities measured in 
arbitrary units have to be corrected by various correction 
factors before being used in calculation of RDF. The major 
corrections arise due to the unpolarized beam used and the 
absorption in the sample. The other corrections are noise 
in the data, multiple scattering in the sample etc. 

The characteristic radiation from X-ray tube is 
considered unpolarized (Bearden, 1928). But such a radia- 
tion after being scattered or diffracted is polarized, the 
amount of polarization depending on the angle through 
which it is scattered or diffracted. The true intensities 
are obtained after dividing by the polarization factor 

P (6 ) = (1 + cos 2©)/2 

The other measure correction is due to the substrate 
scattering. Experimentally recorded intensity Ig c (2©) 
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constitute scattering from both the sample and the substrate. 
To separate out the contribution from substrate, experiment 
is repeated with sample removed from the substrate (inten- 
sity I c (29)). However I (26) does not represent the actual 
contribution from the substrate to I s ^(20) because of the 
absorption in sample. 

Thus the actual experimental intensity from the 
sample is given by 


I expt (29) = I SC (20) “ I c (20) • e" 2 ' MtcOSec9 (2.5) 

where jU is the absorption coefficient of the sample and 
t is the thickness of the sample. 


2.3.3 Intensity measurement 

The sample holder used in the experiment was made 
up of perspex. It was cleaned with acetone and was mounted 
on the stage at the centre of the goniometer. After the 
preliminary adjustments for aligning the diffractometer, 
the step was adjusted at 0.5° in 20. The intensity was 
measured in preset time mode with time preset at 15 seconds. 
Step scanning is activated and intensity in arbitrary units 
(counts per 15 seconds) at each 20 position equally spaced 
in 2© is printed on machine role. The range was from 6° 
to about 130° in 20. The experiment was repeated with 
powder in the sample holder. The tube was operated at 
30 KV, 20 mA. 
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2.4 Experimental Details - X-ray Fluorescence 

When X-rays are absorbed within the sample, a 
bound electron is ejected from the atoms of the consti- 
tuent elements. Electrons from higher levels jump into 
the newly created hole and emit X-ray photon of corres- 
ponding energy. This process is called X-ray fluorescence. 
The fluorescence spectra contains the emission lines of 
the constituent elements, the intensity of the peaks 
being proportional to the concentration. 

The essential features of X-ray fluorescence 
crystal spectrometer are shown in figure 2 ,.4. The 
experimental fluorescence spectra were taken on XRD-6 
spectrometer of General Electric Co. with a fluorescence 
attachment. The primary X-ray beam from X-ray tube with 
tungsten target irradiates the specimen A held in some 
type of drawer B in the specimen compartment C. A lead 
shutter mechanism (not shown in figure) is provided to 
block the primary beam at X-ray tube window when specimen 
is inserted and removed from its compartment. The secon- 
dary radiation thus excited contains many wavelengths and 
is emitted in all directions. But only those rays directed 
substantially parallel to the foils of the solar collimator 
D emerge from specimen compartment and reach the angular 
crystal (LiP)E. Different wavelengths are detected at 
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different 28 angles. The chart recorded v/ith continuous 
scan (1.5° per minute) gives peaks corresponding to the 
emission lines of the elements. The measurements on 

^ e x^ e l-oc s i rs '*' em We £"e used to confirm the concentration 
of constituent elements. 



CHAPTER III 


THEORY AND ANALYSIS OF EXAFS 
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3 . 1 Mechanism of X-ray Absorption 

The experiment of Bragg and Bragg (1913), which 
demonstrated that X-rays could be reflected from crystal 
planes, laid the foundation of X-ray spectroscopy. The 
interaction of X-rays with matter results in a number of 
interesting phenomena, such as scattering, photoelectric 
absorption, radiation damage, pair production (in case of 
X-ray of energy greater than 1.02 MeV), secondary emission. 
Auger emission, etc. The analysis of the attenuation of 
X-ray intensity when X-rays penetrate through matter requires 
the segregation of various physical processes involved of 
which photoelectric absorption is usually the most important 
and is referred to as the true absorption. 

When an x-ray beam passes through a medium the inten- 
sity of transmitted beam I is attenuated logarithmically 
(Azaroff, 1974). According to the classical equation, if 
the incident beam with intensity I q has travelled a distance 
x, then the absorption coefficient will be given by 

= - “ In(f-) (3.1) 

o 

It follows from the above equation that dimensions 
of /t^are reciprocal centimeters so that it represents the 
attenuation of the beam per unit length travelled and is 
called the linear absorption coefficient. It turns out 
to depend on the energy (or wavelength) of the X-rays. 
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Calculation of the linear absorption coefficient is 
complicated by this last fact so that it is convenient to 
define a mass absorption coefficient 


tm = t'/f (3.: 

where f is the density of the absorbing medium,, is the 
cross-section of the mass unit. 


The variation of the mass absorption coefficient 
with the wavelength a. in 51 (10 - ® cm) is shown in Fig. 3.1. 

It is evident from the figure that as the wavelength of 
the transmitted X-rays increases, the photoelectric absorp- 

3 

tion increases as \ until a critical wavelength is reached. 
For wavelengths larger than the critical wavelength at the 
absorption edge, incident X-rays have insufficient energy 
to knock out a particular kind of bound electron in the 
atom, and photoelectric ejection of such electrons no 
longer contributes to the absorption process. The absorption 
edges in Fig. 3.1 are called K, L-j., and edges 

respectively according to the inner electron whose binding 
energy equals the energy at the absorption edge. Following 
the photoelectric ejection of an inner electron, the exeited 
atom typically emits an X-ray photon, as an outer electron 
falls into the newly created hole, and this phenomenon is 
termed X-ray fluorescence. 



■WAVELENGTH (not to seek) — > 

3.1 Schematic L illustration of the variation 
of the aass absorption coef f id ent with wavelength. 
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3 . 2 Extended X-ray Absorption Fine Structure ( EXAFS )_ 

When the horizontal range in Fig. 3.1 is expanded, 
the saw-toothed character of an absorption edge is replaced 
by a fine structure that extends to about hundreds of 
electron volts on the high energy side (or low wavelength) 
of the main edge. This is termed as the Extended X-ray 
Absorption Fine Structure (EXAFS). Absorption process 
refers to the transition of a bound electron to some 
unoccupied final state. The fine structure appearing at 
and near the observed absorption edge (within a few eV) 
therefore reflects the distribution of unoccupied states 
for absorbing atom. As is well Known, proximity of other 
atoms modifies such quantum states and their relative 
densities, so that EXAFS can be utilised in their study. 

This fine stricture was previously observed ( Kronig, 1931) 
only for crystalline substances but later it was observed 
for nearly all types of materials. This was first explained 
by Kronig (1931) for crystalline absorbers. 

3 • 3 Theories of EXAFS 

The first theory explaining the fine structure 
was proposed by Kronig (1931). According to his theory, 
the origin of fine structure is ascribed to the zone 
structure of energy levels of the electrons m the crystal 
which is the periodic field in the crystal. The X-rays 
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are absorbed by electrons which make transitions from X-ray 
levels to empty levels lying above the Fermi surface. These 
empty levels are grouped into zones and there exist forbidden 
and allowed bands. Because the absorption coefficient is 
proportional to the product of transition probability and 
the density of states , a fine structure occurs. In his 
theory variation of transition probability was neglected 
and the fine structure was explained mainly by the electronic 
state density. For the structure of gas molecules he computed 
(Kronig, 1932) the fluctuation of absorption coefficient 
caused by the variation of the state function of the final 
state of the transition. He assumed that final state is 
sum of the plane waves of photoelectron and the scattered 
waves . 

Peterson (1933) applied this theory to explain 
absorption curve of chlorine but subsequent experimental 
measurements by Lindh and Nilson (1943) deviated appreciably 
from his theoretical calculations. Also measurements on 
GeCl^ (Shaw, 1946; Stephenson, 1947; Glaser, 1951) did not 
agree with the predictions of Hartree et al. (1934). 

After failure of Kronig' s theory to explain "general 
features" of the absorption curves of all materials, 

Kostarev (1941, 1949, 1951) presented a more general theory 
based on consideration of nearest neighbour atoms which 
could explain absorption edge structure of GeCl^ (Hartree et al 
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1934). The gist of Kostarev 1 s argument is that the ratio 
■ ¥ - 00 is determined by transition probabilities from initial 
K state to the final P state according to 


T(k) 



2 

free atom 



and is independent of variations in the density of states 
underlying Kronig 1 s theory for solids. 

A different interpretation was offered by Hayasi 

(1949) who suggested that the maxima are superimposed on 

free atom curve and they are additive. Shiraiwa et al. 

(1958) considered the life time of the state of the electron 

ejected by absorption of photon in analysing the original 

Kronig 1 s theory. Assuming that the density of states 
1/2 

varies as E , where E is the energy of photon, they 
concluded that the absorption coefficient 

"X(k) <. e 1/,2 p(e) 


where P(E) is the transition probability. 

It is possible to distinguish these theories from 
early Kronig 1 s theory by realising that Kronig 1 s theory 
(also Hayasi 1 s) is based on the existence of long range 
order (lro) in a crystal whereas other theories consider 
only the nearest neighbour atoms, i.e., the short range 
order (sro). This makes sro theories more general in that 
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they are applicable to amorphous solids and gases as well 
whereas lro exists only in crystals. 

Lytle (1963) began measurements with the initial 
goal of extending Hanawalts temperature effect to cryogenic 
temperatures (Hanawalt, 1931) to sharpen the fine structure 
as well as improving and automating the instrumentation 
(Lytle/ 1966) so that high quality of the data could be 
obtained. Although his theory bears certain similarity to 
the Hayasi theory/ it was successful in extracting accurate 
first neighbour distances from fine structure data. 

Sayers et al. (1970) have reported improved measure- 
ments and theory of EXAFS. Over past few years considerable 
improvement in the theoretical understanding of EXAFS has 
occurred. Theories of EXAFS (Stern, 1974; Ashley and 
Doniach, 1975; Lee and Pendry, 1975) based on the single 
scattering of an ejected photoelectron by atoms in coordina- 
tion shells surrounding the absorbing atom give an expression 
for normalised EXAFS *X_(k) as 

X (k) = ^Ta. (k) sin (2kr j + ct>. (k) ) (3.3) 

j 

where the summation extends over j coordination shells, r. 
is the radial distance from the central absorbing atom to 
atoms in the jth coordination sphere and c£j(k) is a total 
phase shift due to interaction of the photoelectron with 
central atom potential and scattering from atoms in the jth' 
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Aj(k) = (NjA.r^)Pj(k) exp(-2<? 2 k 2 ) (3.5) 

where N. is the number of atoms in the jth shell and <5j is 
the root mean square deviation of the interatomic distance 
about r . which may contain both static and temperature 
effects. The amplitude factor, Fj(k) accounting for elastic 
electron backscattering amplitude f(n ,k) and inelastic 
effect is 

F. (k) = f(rr,k) exp(-2r./k) (3.6) 

being the inelastic mean free path. f ( n ,k) is the EXAFS 
electron backscattering form factor and is a function of 
the type of scattering atom. 

2 2 

— 2 ^ ^ K 

The exponential term e ° accounts for the effect 

of lattice vibration and disorder which was first discussed 
by Schmit (1961) and treated more completely by Beni and 
Platzman (1976). This EXAFS Debye-Waller factor differs 
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from the Debye- Waller factor for X-ray diffraction in an 

interesting wav. It is the mean square average of the 

difference of displacements and may, in general, be 

different from the mean square displacement. The Debye- 

Waller factor has been extracted by studying the tempe- 

2 

rature and the k dependence of the logarithm of the ratio 
of EXAFS amplitudes (Stern et al., 1975), 

The exponential term containing mean free path 
accounts for the inelastic losses in the medium. This 
effect is further dealt with in the following sections. 

Knowledge of backs cattering amplitude and the phase 
shift due to scattering from atoms in the coordination 
shell is essential in comparing the experimental EXAFS 
with expression 3.3. There have been several earlier 
attempts and ab initio calculations of amplitude and phase 
shift by Hartree-Fock (HF) method (Lee and Pendry, 1975; 
Kincaid and Eisenberger, 1975; Pettifer and McMilan, 1977) 
and by Hartree-Fock-Slater (HFS) (Ashley and Doniach, 1975; 
Lagarde, 1976). Neither scheme provided very satisfactory 
results. What is needed is a relatively simple scheme 
which intsrpolates between the low energy region where HFS 
method has been tested, and the high energy region where 
HF approach is expected to be more reasonable. Such a 
scheme based on the local density functional formalism 
was introduced by Lee and Beni (1977). Teo and Lee (1979) 
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have calculated amplitude and phase functions for nearly 
half of the elements in the periodic table using an 
electron-atom scattering theory of Lee and Beni (1977). 
These values have been used in the analysis and found 
to give satisfactory results (Shulman, et al. 1978; 

Teo et al./ 1977; Lee et al., 1911 , Teo et al. , 1978a/ 
1978b). 

3.4 BXAFS Analysis 

One of the earliest attempts of successful extrac- 
tion of structural parameters from the EXAFS data was based 
on the Fourier transform of the data expressed in momentum 
space (Sayers et al., 1971). The absolute value of the 

transform was found to peak at distances shifted from known 

o 

values by several tenths of A . When corrected for these 

shifts using systems with known distances/ bond lengths 

could be extracted. However, this method is not very 

accurate because the peak value of the Fourier transform 

depends in general on such factors as the weighting of the 

data and the choice of the threshold energy E and because 

o 

the phase shift k ) , which is largely responsible for the 

shifted Fourier peak, is not a linear function of k. 

Citrin et al. , (1976) extracted the phase function (k) 

by direct curve fitting to eqn. 3.3 to show that these 
functions could be transferred to other simple systems to 
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o o 

extract bond lengths with an accuracy 0.01 A to 0.02 A. 

This approach, however, is not generally applicable to 
systems of more than one principle distance. The two 
methods of Fourier transform and curve fitting have been 
combined to analyze more complicated systems. Following 
sections provide a relatively detailed account of the 
Fourier transform filtering method which has been used 
in the present investigation. 

3.4.1 Extraction and normalization of EXAFS 

The removal of oscillatory part (EXAFS) of the X-ray 
absorption coefficient from the smooth mo no tonic ally decrea- 
sing background is the most difficult problem as the smooth 
absorption of an isolated atom ^ (k) is not generally 

available experimentally, and the present theoretical 
calculations of /^(k) are not sufficiently accurate for 
most EXAFS work, (nearly .01%). There are several background 
removal methods that could be tried such as simple polynomial 
fit over the whole range of data, iterative low order poly- 
nomial fitting combined with fitting EXAFS, orthogonal poly- 
nomial fitting, Fourier transform, extrapolation of an assumed 
functional form of the background from below the absorption 
edge and the most flexible cubic spline fitting. In our 
analysis the EXAFS signal was separated from low frequency 
background and high frequency noise by Fourier filtering 
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technique (Lytle et al. 1975). A sixth order polynomial 
fitted to the data (EXAFS + background) over the whole 
range approximated /a (K) and used to obtain the normalized 
EXAFS . 

3.4.2 Fourier transform filtering 

The next step in the analysis is a Fourier transform#. 

To perform the Fourier transform, the data must be transformed 
from the energy to wave number k-space. The transformation 
is given by 

^ 2 k 2 /2m = E-E q (3.7) 

where E is the energy of the X-ray photon and E o is the 
"inner potential" caused by the atomic potentials and 
represents the "zero" above which kinetic energy must be 
added to determine the total energy E. In metals at low 
energy, E q can be approximated by Fermi energy, but at higher 
energies ( > 30 eV or so) this is no longer true^ since the 
exchange and correlation energies go to zero in a complicated 
way with increasing energy, and the pseudopotenii al of the ion 
core is energy dependent. 

Several workers have fixed E q arbitrarily at a special 
feature in the spectra such as a point of inflexion 
(Hayes et al.,1976 ; Cramer et al., 1978). The effect of 
the choice of E q and its determination are discussed in 
the later part of the chapter. 
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To make the Fourier transform less sensitive to the 

3 

choice of E , the normalized EXAFS is multiplied by k 
o 

3 

and then Fourier transformed. The k multiplication has 
the effect of weighting the EXAFS oscillations more uniformly 
over the range of data (k >4A “ 1 ) and also assures that the 
chemical effects on EXAFS information, which are most 
significant at small k are effectively cut off. The magni- 
tude of Fourier transform peaks at several distances shifted 
towards origin by an average slope of phase shift <p (k) . 

The first shell contribution to EXAFS can be isolated 
by applying a smooth filter window function. Both the real 
and imaginary parts of Fourier spectrum are filtered identi- 
cally so that the phase error is minimised. The essential 
feature of a good filter window function is that it does 
not produce gross distortions in the data and almost any 
smoothly varying function satisfies this requirement. The 
first shell EXAFS can be decomposed into amplitude and 
phase functions using the Fourier transform method 
(Shulman et al., 1978). 

3.4.3 Empirical phase shifts and chemical transferability 

The argument of sine function in eqn. 3. 3 is 
-^p.(k) =2 kr . + ( p ■ ( k) which can be obtained experimentally 
for a model system. At sufficiently high photoelectron 
kinetic energies, typically above 50 eV where EXAFS 
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scattering processes are dominated by core electrons/ phase 
shifts are insensitive to chemical environment. This is 
the principle underlying the concept of chemical transfer- 
ability in which it is assumed that the phase shift in the 
case of model system is the same as that in the case of the 
unknown system. Determination of'vpj(k) for an atom pair 
in a given model system of known distance rd enables the 
determination of distance r^ in an "unknown" (different) 
system containing the same atom pair as is evident from 
the following. 

For the model system/ the phase function is 
^(k) = 2kr m + 4>(k) 

where ‘dp (k) is the scattering phase shift defined earlier. 

For the unknown system containing the same atom 

pair 

^d(k' ) = 2k'r u + C|>(k) 

using the concept of transferability of phase shift. 

The wave vector k' for unknown system is not/, in general/ 
same as k of the model due to the chemical shifts/ i.e. , 
because of the change in the threshold E q . 

The plot of -^(k) - ^(k 1 ) as a function of k 

should be a straight line with zero intercept and slope 
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equal to 2 (r m -r u ) . E is varied in such a fashion so as 

o 

to reduce the intercept of the line to minimum. A closely 
related analysis method of varying "E until the "slope" 
remains constant over a large range of k has been used 
successfully by Martens et al. (1978). 

Another flexible method which has been used in our 
analysis, of determining E q and r simultaneously, is the 
least square fitting the experimental phase shift 4>oo 
with the theoretical phase shift (Teo and Lee, 1979). Since 
tj>(k) = ^ (k) - 2kr, 

£s- r 4 > = 2 k <Lr 

Thus changes in r affect (k) mainly at higher k values. 
Also A. *3^ = 0.2 63 AE /k indicates that changes in E 

affects only at lower k values. Thus it is not possible 
to produce an artificially good fit by adjusting E q with 
a wrong distance r. This method has been widely used and 
the nearest neighbour distances could be determined to 
within 0.03 % (Reed et al., 1977, 1978; Teo et al., 1978a, 
197 8b) . 

3.4.4 Amplitude transferability 

The amplitude function obtained after decomposition 
of first shell EXAPS contains information about the coordi- 
nation number, the kind of atoms involved, the amount of 
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disorder and the effective photoelectron mean free path, 
In principle/ it is possible to obtain all these para- 
meters by the direct method of curve fitting. But in 
practice, the experimental amplitude is usually off by 
50%. This discrepancy is expected to be somewhat 
sensitive to the chemical environment and will depend 


on distance r. (due to an exponential factor e 


-2r j/X 


Moreover, reduction in amplitude is due to the many body 
effects (Stern et al., 1980) which are nob considered in 
the theory. The EXAFS amplitude is also affected by thick- 
ness and non-uniformity of the absorber (Stern and Kim, 1981), 
Because of these reasons and the strong correlation between 
the various parameters determining the amplitude function, 
the curve fitting technique is not suitable. 


Analogous to the chemical transferability of phase 
shifts, in principle, it is possible to extract these 
quantities assuming chemical transferability of amplitudes, 
i.e., using model compounds as standard for unknown systems. 
For a single shell containing N identical atoms, the 
amplitude function is given by 

A(k) - " ,(k> e- 25 - 2 * 2 e- 2r ^ (k > (3.8) 

kr z 

Denoting the model and unknown systems by m and u 
respectively, the logarithm of their ratio gives 
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If it is assumed that X (k) 25 -^* (k) and (r -r ) ^V'XCk), 

then the third term becomes negligible and the parameters 

of the straight line obtained by plotting the logarithmic 

o 

ratios of the amplitude versus k enable, the evaluation of 
and . A similar analysis could be used for deter- 
mining N. for different shells. 

There are several difficulties that can be encountered 
with this approach. The first involves empirical determina- 
tion of A(k). The finite value of k (the maximum value 

•max 

of k) affects the magnitude of the Fourier transform and in 
turn affects the amplitude function than its phase. Another 
factor is the relative mean free path of the two systems. 
Differences as little as 1-2 A° at a given k are actually 
non-negligible because of the logarithmic form of the above 
equation (citrin et al. # 1979). This problem has also been 
discussed by Eisenberger and Lengeler (1980) and Stern et al. 
(1980) . 



CHAPTER IV 


CRYSTALLINE AND AMORPHOUS GERMANIUM 
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4.1 Introduction 

* 

In this chapter are presented the results on 
crystalline and amorphous germanium. Both crystalline 
and amorphous germanium have been studied before with 
X-ray diffraction (Shevchik and Paul, 1972; Temkin et al., 
197?) and EXAF3 (Sayers et al. , 1971, 1972? Crozier and 
Seary, 1981). Earlier work of Sayers et al. (1972) demons- 
trated that EXAFS could provide quantitative information 
regarding the local structure around the absorbing atom. 

In view of the progress made in theoretical aspects and 
data analysis procedures in the last decade, it is worth- 
while to study crystalline and amorphous germanium with 
EXAFS. Experimental data on crystalline germanium (c-Ge), 
which serves as model system in the analysis of amorphous 
germanium (a-Ge), has been used to standardize the analysis 
procedure. Although information about second and higher 
coordination shells is important in distinguishing different 
structural models for disordered systems, first shell para- 
meters provide a useful aid. The different data analysis 
procedure for ne-r neighbour distance determination, in 
particular, have been examined. 

4.2 Elementary Analysis 

The absorption j*K x as a function of incident X-ray 
phton energy for K-absorption discontinuity in c-Ge is 
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depicted in figure 4.1. The abrupt change in the absorption 
near the absorption edge around 11 KeV amounts to nearly 1.5. 

A sharp peak in the vicinity of edge on the high energy side 
is the white line (WL). It may be compared with the earlier 
work of Lytle et al. (1975). The broadening of the white 
line observed may be attributed to the scattering effects 
and relatively higher thickness of the absorber (Agrawal, 
1979). The broad peaks after the WL vanishing far above 
the edge are due to EXAFS which are clearly visible in the 
figure. The slowly decreasing absorption on the low energy 
side of the absorption edge is due to the substrate and any 
other contribution resulting from the other bound electrons 
in the atom. This was approximated by a Victoreen curve of 
the form A X 3 + bA 4 with A = 0.349; B = 0.012 4 (shown by 
dashed curve in figure 4.2(a)). The contribution to absorp- 
tion exclusively due to K shell was obtained by extrapolating 
and subtracting this curve from the data on the high energy 
side. 

The wave vector at every point above the edge X > 60 eV) 
was calculated using the equation 3.4. The value of energy 
at which the kinetic energy of the photoelectron is zero 
was obtained by comparing the slopes of the absorption edge 
at every point on it. The average of the energy of the point 
at which the slope starts decreasing and that of the prece- 
eding point was chosen as the zero point energy (E = o). 
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This approximately corresponds to the inflection point of 

the edge (11129 eV). Since the wave vector is non linearly 

connected with the Bragg angle 9, the k values are not, 

equally spaced. A quadratic interpolation programme (see 

Appendix I ) has been used to interpolate the data onto a 

. o_i 

grid of equally spaced k values (h>k~0.02 A ). The 
resulting data now serves as input to the Fourier transform 
routine in the next step of analysis. 


4.3 Removal of Free Atom Absorption fr xj. k)and Normalization 


One of the most difficult problems in the EXAFS analysis 

is the removal of background and the normalization of EXAFS 

oscillations to per atom scale. Direct Fourier transform 

of the raw data after adding zeroes from 0 to k . and from 

a mm 

k m in to about 12 k indicated that the background frequencies, 
(low frequencies) extend to the intermediate EXAFS frequencies. 
The EXAFS frequencies could be easily estimated because the 
first shell radius for crystalline germanium is known to be 
equal to 2.45 %. Therefore all the Fourier coefficient 
(a^s and b^s) corresponding to r values greater than ^1.3 %. 
to about 7 A represent the EXAFS signal. A sixth order 
polynomial in wave vector k, which also approximates (k) 
was fitted to the data (EXAFS + background). Fourier trans- 
form after separating /-*• (k) gives EXAFS frequencies well 
separated from the low frequency background. EXAFS in 
k-space was obtained by performing an inverse Fourier 
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transform (IPT) after removing low frequencies corresponding 

to any residual background and the high frequencies corres- 

o 

ponding to r greater than 7 A representing noise. This was 
then' normalized to per atom scale using the polynomial fit 
( (k ) ) . 


Figure 4.2(b) shows the resulting normalized EXAFS 

X(k). As is evident from the figure, EXAFS peaks around 
°_1 

3 A and decays at high k-values almost exponentially. 
Relatively high value at low k is due to the distortions 
produced in the process of extraction and normalization 
of EXAFS and such distortions are not desirable for further 
analysis. The difficulty is resolved by multiplying yj> k) 
with a Hanning window function (Bergland, 1969) defined as 


\ \ 1 --cos2^(k- k . )/(k 

2 ! ' mm max 


mm 




This function has been used in the region of k-values 

corresponding to the first and last 10% of the range investi 

gated. As a result 'V(k) is forced to assume a half cosine 

bell shape terminating at k . and k _ (Lytle et al. 1980) 

mm max J 


4. 4 Fourier Transform Filtering 

The normalized EXAFS thus extracted is, in principle, 

capable of giving the local structure parameters such as 

radii of different coordination shells, corresponding coor- 

2 

dination numbers and the disorder parameters C . etc. by the 
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method of least squares fit using equation 3.3. However 
the list of parameters is so large that such a least squares 
minimization is practically not possible. Another flexible 
method briefly described in the previous chapter, namely 
the Fourier transform filtering, has been employed here. 

Figure 4.2(c) shows the magnitude of Fourier transform 
(FT) of k "X(k). The k weighting of the SXAFS data approxi- 
mately equalizes the amplitude over the whole range and 
weighs the data at large k values properly. This further 
makes the position of the peaks in the FT less sensitive 
to the choice of E . Recently it has been shown (Stearns, 

19 82 ) that at a particular choice of E q the phase shift 
■k 

Cp(k) is linear and the position of peak in the FT is 

3 5 

independent of the weighting factor (k, k , k ). This will 
be discussed later. 

o 

The first large peak is observed at 2.27 A and 

corresponds to the contribution of the first coordination 

shell to EXAFS. The other prominent peaks correspond to 

o 

the higher coordination shells. Structure beyond 6 A is 
not visible as it was filtered out during the extraction 
of EXAFS. The dashed curve encompassing the first peak 
is the smooth filtering window which was used to separate 
out the first shell EXAFS. The resulting first shell 
EXAFS, k 'X^(k) obtained after performing I FT is shown 
in figure 4.2(d). If before performing IFT, the coefficients 
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corresponding to the negative r— values were removed, the 
resulting inverse transform would be a complex quantity 
(see Appendix n ) from which the amplitude and phase 
functions could be readily evaluated. The amplitude 
function k 'A^(k) obtained in such a way is shown by the 
dashed envelope in figure 4.2(d). The phase function 
(k) = 2kr^ +<d^(k) is shown in figure 4.3. As <j>^(k) 
is very small compared to the factor 2kr^, the curve 'vL^ (k) 
is almost a straight line with slope equal to 2 (r^ - c<) 
where ©c is called the bonding parameter ( 4* (k) = ~2»kk-;-jS ). 
The peak in the FT is shifted by an amount ©c towards the 

origin. The least squares straight line fit gives value 

o o 

of = 0.19 A using the known value of r^ = 2.45 A. This 

may be compared with the earlier work of Lytle et al.(l975) 

o 

where the value of c J. is equal to 0.2 4 A. The discrepancy 

is due to the value assigned to E . Change in value of E q 

changes the position of the peak in FT and subsequently 

the value of &<. obtained. Thus any physical significance 

to ©< could be assigned only when certain criterion for 

the choice of E is fixed. 

o 

4.5 Linear Phase Shifts and Bond Length Determination 

As noted earlier the weighing of the EXAFS data 
reduces the sensitivity of position of the peak in the 
magnitude of FT on the choice of E . 


The fact that the 
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peak position depends on the weighting factor as well as 

on the choice of E q has been posing a problem in determining 

the accuracy with which bond lengths could be determined 

with EXAFS. Recently (Steams, 1982) it has been pointed 

out that there exists for any material a particular value 

of E q = E^ at which the peak position in the magnitude of 

FT is independent of the weighting factor (k, k 4 k ) and 

that the phase shift derived at this value of E q is linear. 

For a particular weighting factor, the variation of peak 

position as a function of E q is nearly a straight line, 

the slope of which diminishes with increase in degree n of 

the weighting factor. The lines for different values of n 

intersect at a particular point which gives the value of 

E . Further the existence of such a crossover point can 
c 

be used to assess the quality of the experimental data. 

Using the method outlined above, the value of E c 

has been determined for c-Ge as well as for a— Ge (see 

section 4.7). Figure 4.4(b) shows the variation of the 

position of first peak as a function of E q and weighting 

factors k, k 3 and k 5 for c-Ge. The crossover point obtained 

crives value of E = -2 5 eV and is comparable with E = -2 6 eV 

c u 

as obtained by Stearns (1982). The near neighbour distance 
can now be determined from the phase function obtained at 
E = E c by comparison of experimental and theoretical phase 
shifts. Figure 4,5 shows the two phase shifts obtained 
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after the minimization procedure. The dashed curve is 
the theoretical phase shift (Teo and Lee, 1979). The 

o 

parameters determined are AE = 30 eV, r 1 = 2.46 i 0.01 A 

° o 

(with width of the filtering window equal to 1.2 5 A). 

4.6 Effect of the Width of the Filtering Window 

The effect of the width of the filtering window on 
the phase and amplitude function and subsequently on the 
determination of local structure parameters has been 
investigated by analysing the normalized EXAFS data with 
different widths of the filtering window. The phase func- 
tions were used to determine the near neighbour distance 
by comparison of the experimental and theoretical phase 
shifts. The changes in the phase functions could be 
immediately ascertained from the values of and the near 
neighbour distance obtained for different widths of the 
filtering window. The results are summarised in Table 4.1. 

The value of r 1 determined by comparison of phase shifts 

o 

with different widths is same and equals 2.46 A with the 
accuracy of ± 0.01 A. The value of «< or slope of the phase 
function does change with the width of the window. However 
this does not affect the value of r^. 

Unlike the case of phase function, the width of the 
filtering window not only affected the magnitude but also 
the shape of the amplitude function. This is illustrated 
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Table 4.1 


Effect 

of width 

of the filtering window 

on phase function 



(E o = 

-2 5 eV) 


Width 

(A) 

(%) 

A.E 0 (eV) 

r ± d) 

0.535 


-0.058 

34.05 

2.46 

0.714 


-0.059 

34. 04 


0.892 


-0.060 

33.80 

2.46 

1.071 


-0.058 

33.11 

2.46 

1.250 


-0.042 

30.20 

2.46 

1.428 


-0.013 

25.00 

2.47 


1.607 


-0.015 


19.70 


2.47 
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by plotting the amplitude function A^ (k) for five diffe- 
rent widths of the filtering window in figure 4.6. As 
is evident from the figure a too narrow window has comple- 
tely washed away the details in the amplitude function. 

With increasing width all the amplitude functions start 

°-l 

with a lower value and peak around 7 A (except for 
'narrow' window). Thus for sufficiently large window 
width it is only the low k information which is mostly 
affected. The theoretical calculations of backscattering 
amplitudes (Teo and Lee, 1979) suggest that the amplitude 
function for Ge has a peak around 7-8 i 1 and then decays 
at large k values. For extracting the structural informa- 
tion from the experimental amplitude by comparing it with 
the theoretical one (equation 3.8’) , the two functions must 
be consistent. Non-uniqu--ness of the amplitude function 
naturally causes concern and great caution should be 
exercised while determining the structural parameters from 
this function. In our analysis we have used an optimum 
width (1.25 A) for which the details in the amplitude 
function are retained. 

4 • 7 Results for Amorphous Germanium 

Similar analysis has been carried out for the "K 
absorption discontinuity of Ge in amorphous germanium 
( a _Ge). The results are presented in graphical form in 
figure 4,7. Figure 4.7(a) shows the experimental absorption 
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spectrum with dashed curve representing the Victoreen fit 
(A = 0.9319; B = -0.1749),. The steep decrease in absorp- 
tion as compared with corresponding spectrum for c-Ge is 
duo to the substantial contribution to absorption by the 
substrate (absorber was prepared by sandwitching several 
thin films together). The normalized EXAPS obtained from 
raw data after background removal and subsequent normaliza- 
tion to per atom scale appears in part (b) of the figure. 

The nature of EXAPS oscillations is similar to that obtained 
in the case of c-Ge except the rapid changes in the negative 
half. The reason for this could be explained by considering 
the absorption spectra of c-Ge and a-Ge. The absorption 
spectra of a-Ge contains a lot of noise as compared with 
that of c-Ge spectrum. Since the noise frequencies are 
uniformly distributed throughout the entire frequency range, 
it is not possible to separate the noise completely. 

3 

Figure 4.7(c) shows the magnitude of FT of k '^(k) 

with dashed curve being the 1 filtering window. The first 

o 

large peak is observed at 2.31 A. The overall spectrum 
is similar to the one obtained for c-Ge. In figure 4.7(d) 

3 

is depicted the first shell EXAPS k ‘')L(k). It rises 

°-l 

sharply and peaks around 12.5 A . The dashed envelope 

3 

is the amplitude function k A^(k). The phase function 

( k) gives the value of <*: (calculated using r^ = 2.45 A) 

o o 

equal to 0.16 A, lower than that obtained for c-Ge (0.19 A). 
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This deviation of nearly 0.03 A in the value of °< indicates 
a larger value of r^ in a-Ge. The position of the first 
peak in the magnitude of FT which is shifted by an amount 
towards origin gives r 1 =2.$0-S. for a-Ge (using o4 of c-Ge). 

The value of E c for a-Ge was determined by the proce- 
dure described in section 4.5 and turns out to be -55 eV 
(fig. 4.4(a)). The first shell radius now can be determined 
from the phase function extracted at E c = -55 eV. Figure 4.8 
shows the comparison of experimental and theoretical phase 
shifts obtained after minimization. The parameters obtained 
are r 1= 2.49 A, ^E q = 59 eV. The value of r 1 obtained by 
this method agrees well with the straightforward method of 
usings for c-Ge (at E q = 0). 

4 * 8 Empirical Phase Shifts and Chemical Transferability 

The concept of chemical transferability has already 
been outlined in section 3.4.3. Accordingly a difference 
in phase functions of c-Ge (model) and a-Ge was calculated. 

A least squares straight line fitted to this had intercept 
-2.31. Non-zero value of intercept indicates that the 
threshold in a-Ge and c-Ge differ. Z1 E q was varied using 
the method of bisection till the intercept reduced to 
nearly zero. Using the predetermined value of for 
c-Ge (2.46 ± 0.01 A), value of slope of the straight line 
yields r ± = 2.49 ± 0.01 A for a-Ge. Figure 4.9 shows the 
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% 

plot of the difference of phase functions and the least 
squares fit. The straight line clearly passes through 
origin. 

4 * 9 Determination of Coordination Number 

Information about the coordination number is con- 
tained in the amplitude function obtained by FT filtering. 
Comparison made between measured EXAFS amplitudes for 
nearest and more distant neighbour coordination shells in 
wide variety of known compounds at room temperature and 
liquid nitrogen temperature reveal the significant sensi- 
tivity of amplitudes to chemical environment (Eisenberger 
and Lengler, 1980). Absence of exact quantitative theore- 
tical calculation of amplitude function and systematic 
errors in the experimental data makes the determination 
of coordination number by direct comparison of amplitude 
function unreliable. However an overall sealing factor 
determined for model system may be used, to a first 
approximation, in the other chemically, similar system to 
determine the coordination number. The other widely used 
method is chemical transferability of amplitudes with known 
chemically similar model system (see section 3.4.4). 

4 . 9 . 1 Comparison of amplitude functions 

The amplitude function for the first coordination 
shell is given by expression 3.8. Mean free path A can 
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be assumed to be independent of the photoelectron wave 
vector k. Using the coordination number for c-Ge system 
(4.0) and the quantitative theoretical values of back- 


scattering amplitudes (Teo and Lee, 1979), an overall 

2 

sealing factor and Debye waller type f actor ^ could be 
determined from the experimental amplitude function for 
c-Ge. Figure 4.10 shows the final result of the minimiza- 
tion run. Dashed curve is the theoretically calculated 
kA 1 (k) curve. The parameters namely the overall scaling 
factor and cT 2 are 0.2336 and 0.0034 A* 2 respectively. 

Using these parameters in fitting the amplitude function 
for a-Ge with theory, the coordination number for a-Ge was 
determined and is equal to 4.4 + 0.8. In figure 4.11 the 
two functions kA.^k) for the case of a-Ge have been 
compared. 

4.9.2 Amplitude function s and chemical transferability 

Assuming that the two amplitude functions for c— Ge 
and a-Ge are transferable (i.e. backscattering amplitudes 
F(k) remain same in both), logarithmic ratio of the two 
functions was plotted against the square of the wave vector 
k. However the plot is no more a straight line as could 
be expected from the expression 3.9 (after neglecting the 
third term) . Therefore a least squares fit with a straight 
line was obtained. It was found that the slope and the 
intercept of this line vary depending on the range of 
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k— values over which the fit is obtained* Therefore the 
parameters namely the coordination number for a-Ge and * 
tne difference in ^ values for c-Ge and a-Ge could be 
determined only approximately. Figure 4.12 shows the 
logarithmic ratio and the fitted straight line plotted 
against k over a range 5 — 15 A t The slope' and 
intercept thus determined are -0.00048 and -0.072 respec- 
tively. Neglecting the third term in equation 3.9., the 
coordination number for a-Ge turns out to be 4.49 ± 0,5. 

4*10 Results : X-ray Diffraction 

X-ray diffraction pattern of a-Ge films show 
broad halos and absence of any sharp peaks which indicate 
the amorphous nature of the films. The e:xperimental data 
consisted of measured intensity from the Al-substrate and 
the film (I q and I respectively). The contribution to 
intensity I from the Al-substrate was subtracted using 
equation 2.5. The intensity in arbitrary units was further 
corrected for polarization, compton scattering and dispersion 
of form factors. The normalized intensity is obtained by 
requiring that intensity curve at high s (4 Asin0/ 
oscillates about the total independent scattering curve 
(Haider et al„, 1966). Figure 4.13 shows the normalized 
intensity curve. The normalized intensity I au is related 
to the radial distribution function as follows : 
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F(s) = 


I 

( eu „ ] 

si ?~ 1 = 


CP 


i 


G(r) sin sr dr 


where f is the atomic scattering factor, r is the radial 
distance from an atom at an arbitrary origin; G(r) equals 
4?cr( f - f Q ) and is the Differential Radial Distribution 
Function (DRDF); s equals 4 7fsinQ/x ; f is the atomic 
density, atoms per unit volume, at a radial distance r, 
and ^ is the average atomic density. 


Fourier inversion of the above equation yields 


G(r) = ~ j F(s) sin rs ds 

Following the termination correction procedure 

(Kaplow et al. , 1965), the final G(r) function was obtained 

and is shown in figure 4.14. The near neighbour distance, 

o 

from the first large peak, is 2.46 ± 0.03 A. This is in 
agreement with the value obtained from EXAFS analysis. 



CHAPTER V 


Se„ SYSTEM 
X 1-x 
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5.1 Introduction 

Ge^Se^_ x compounds can be obtained in crystalline 

and amorphous structure depending on the concentration 

of germanium (x = 0.33 and 0.5, the compounds have 

crystalline structures). The system Ge Se, system has 

xl -x x 

been studied earlier by X-ray diffraction methods 
(Fawcett et al., 1972? Bienenstock,i97'^b) and also by 
EXAFS (Sayers et al. , 1972, 1974) and Raman spectroscopy 
(Nermanich et al., 1978; Kawamura et al.,1982). The X-ray 
diffraction studies, in particular, met with several 
limitations. For low concentration of germanium (x^-0.33) 
both the chain crossing model (CCM) and Random covalent 
model (RCM) were equally probable. Earlier studies with 
EXAFS on amorphous GeSe and GeSe 2 too could not point at 
a particular model which would describe the experimental 
results completely. The present work on this system is 
limited to compounds with low germanium content. 

In this chapter are presented the results obtained 
for two compositions x = 0.1 and x = 0. 3 (hereafter 
referred to as I and II respectively). The experimental 
data of both X-ray diffraction and X-ray absorption (EXAFS) 
experiments have been analysed with more emphasis on later. 

The K X-ray absorption discontinuities of germanium 
as well as selenium have been recorded in both I and II. 
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Analysis of EXAPS associated with particular kind of 
absorbing atom provides information about local structure 
around that atom. However the proximity of atomic numbers 
of Ge and Se poses some problem in determination of near 
neighbour distance using experimental phase function. An 
entirely different approach has been used/ for the above 
mentioned reason, which leads to a model based on probabi- 
lities of all the possible configurations. The experimental 
first shell EXAFS could then be described in terms of this 
model. 

5.2 Results ; EXAFS 

In all the cases the elementary analysis and extrac- 
tion and normalization of EXAFS has been carried out by 
procedures described in the previous chapter. The results 
for Ge K absorption discontinuity in I and II are presented 
graphically in figures 5.1 and 5.2 respectively. The fine 
structure is clearly visible in both the absorption spectra 
(part (a) of the figures). Incident photon energy at the 
inflection point (E =0) in both cases is 11126 eV. 

The background absorption below the edge is more in case 
of I. This is because the sairple of same thickness could 
not be prepared with the particular method used for pre— 
paring the absorption screens (see section 2.2.5). Tne 
normalised EXAFS shown in (b) part of the figures is 
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comparatively sharply peaked for the case II. However 

the overall trend in both the cases is similar starting 

with low amplitude and giving a peak around 7 The 

magnitude of Fourier transform of k^”)C(k) (with E q = 0) 

appears in figure 5.1(c) and 5.2(c). Apart from a sharp 

main peak corresponding to first coordination shell, a 

not well resolved small peak at 2.6 A is observed in case 

of I. Such type of: structure has been previously observed 

by Sayers et al. (1972) while studying a-GeSe, However 

this type of structure is absent in II. The position of 

the main peak in II is at slightly higher distance 
. o 

(“'-' 0.03 A) than that in I. The other striking feature 
is the presence of sharp peak at 3.37 A in I which is not 

clearly resolved in case of II. The first shell EXAFS 

3 A 

and the amplitude function k ^(k) are shown in 

figures 5.1(d) and 5.2(d) for I and II respectively. 

o 

With ^the same width (1.25 A) of the filtering window used 

for filtering out the first shell EXAFS, the amplitude 

°-l 

function peaks around 10 A 

Similar results of EXAFS associated with K X-ray 
absorption discontinuity of selenium in I and II are shown 
in figures 5.3 and 5.4 respectively. The absorption spectra 
shown in part (a) of the figures reveal that the background 
absorption due to the sample is low in the energy range 
12-14 KqV. The inflection points are observed at 12776 eV, 
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127 71 eV in I and II respectively. Unlike the case of Ge 
K edge, the EXAFS oscillations are not prominent in this 
case and also the white line is absent. As is evident 
from the figures the patterns of normalised EXAPS oscil- 
lations are also very different when compared with the 

case of Ge K discontinuity. The first large peak for I 
o 

(2.12 A) is at slightly higher distance compared to the 

o 

position of the corresponding peak in II (2.03 A). As 
in case of Ge K edge, the splitting of the first main 
peak is observed in this case also both for I and II. 

This splitting which is more clearly resolved around Se 
atoms may illustrate the tendency of Se atoms to link two 
nearby Ge atoms into a chain as the other Ge atom is forced 
away. 

Although the structure observed beyond first peak 
could be of considerable significance, here the analysis 
is extended only up to the first coordination shell. 

5.3 Determination of Near Neighbour Distances 

The determination of near neighbour distance from 

experimental phase function by comparing theoretical and 

experimental phase shifts reguires that the phase function 

be extracted using E = E . The critical value E for 

o c c 

all the cases (four) was determined by the procedure 
described in section 4.5. Table 5.1 summarizes the 
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Table 5.1 


s 

c 

values in I 

and II 

for first peak 


Compound 

Absorbing 

atom 

E c 

(eV) 

Peak position 
at E c (&) 

Peak position 
at s o o= o; ( i) 

I 

Ge 

-50 

2.230 

2.040 

I 

Se 

5 

2.100 

2.120 

II 

Ge 

-10 

2.050 

2.017 


II 


Se 


42 


1.820 


2.030 
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results of E^ determination. Moreover the calculation 
of theoretical phase shift which is given by 

( k) — ^ — ~JT 

where c£> a is central atom phase shift and is the 

backscattering atom phase shift, requires that the atom 
pair i.e., the central or the absorbing atom and the back- 
scattering atom is identified. 

Unlike the case of crystalline and amorphous 
germanium where the backscattering atoms in the coordina- 
tion shell are of same kind, in case of Ge Se- system 

X -L —X 

due to the proximity of atomic numbers of Ge and Se, either 
Go or Se can act as the backscattering atom. - For example 
consider the first shell phase function obtained from EXAFS 
associated "with Ge K absorption discontinuity in I. The 

•ti*’ 

experimental phase function for the first shell represents 
the resultant of phase function of possible bonding of 
central Ge atom with another Ge atom and Se atom. However 
the two possible bondings may not occur with equal probabi- 
lities cs the concentrations of Ge and Se are very different. 

One can assume that the backscattering atom is a 
germanium atom and determine the near neighbour distance 
r-L and change in threshold A E q by fitting the experimental 
phase shift with theoretical one. Figure 5.5 shows the 



m 



QsmpmriMon of tu^arlinnntftl mi 6 th«o- 
£wx.le*l (dach«dj pfra«« shifts of 0* X mdgm In 
I* BackJDoat taring atom is 



99 


two phase shifts obtained after the minimization run. The 
dashed curve is the theoretical phase shift. The sum of 
squares of the deviations obtained after this fitting is 
a measure of the quality of fit. On the other hand the 
fit would be better if Ge itself (assumed to be backs cat tering 
atom ) is the more probable neighbour. Denoting the sum 
of squares of deviations by (F) , where x is the kind of 

X 

backs catte ring atom, the probability P of x being the 
nearest neighbour is inversely proportional to (F) x * This 
probability would also depend on the concentration of 
the X-atoms in the compound. 


Thus we have 


or 


p x 


I< C 


TfJ 


X 


whore K is the constant of proportionality which in general 
will depend upon the quality of experimental data itself 
but not on the kind of backscattering atom. 


For the present case of Ge edge in I 


K C 


'Ge 


0.1 K 


'Ge 


(P) 


Ge 


(F) 


Ge 


(5.1) 


Similarly one can assume that the backscattering 


atom is Se and obtain 
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K C, 


Se 


Se 


(F) 


0.9 K 

TfT 


(5.2) 


'Se Se 

Here K remains same as the same experimental phase 
function is under consideration. 

Since the bachscattering atom is either Ge or Se, 
one must have 


P Ge + P Se “ 1,0 


(5.3) 


Using equations (5.1), (5.2)/ (5.3) the probabilities 

p , and P n could be determined. Around Ge and Se in both 
Ge Se 

I and II we have in all eight different cases to be consi- 
dered. Similar calculations have been carried out for all 
those cases and results are tabulated in Table 5.2. The 
probability P x appear in the last column of the table. The 
vr.lU'S oE Ae obtained are consistent with respective E c 


vc luos . The near neighbour distance determined for I 

corresponds to covalent bond between Ge and Se. However 

o 

for II the value of r^ is smaller by — 0.04 A. 


From the probabilities for all the cases it follows 
that Se is th- most probable nearest neighbour around Ge 
W j 11 as Se. Also there is a finite probability of Ge 
boil g nearest neighbour of Ge and this casts some doubt 
on tho chain crossing model proposed earlier (Fawcett et al 
1972). In the so called "chain crossing model (CCM) ,f , the 
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Table 5.2 

Probabilities P and P c obtained 

Cos be 

from phase functions 


Compound Central 

1 atom 

Sackscat- 
tering atom 
_ (x) 

r l (A) AE (eV) 
o 

(3?) x 

p. c CO 

Ge 

I 

Ge 

2. 378 

57.49 

0. 0121 

13.9 

Ge 

Se 

2.370 

61.33 

0. 0176 

86.1 



Ge 

Ge 

2.338 

7.23 

0.0415 

27.0 

ri 








Ge 

Se 

2.328 

10.88 

0.0370 

73.0 


Se 

Ge 

2.387 

12.48 

0.0293 

8.5 

Se 

Se 

2.37 8 

16.00 

0. 0249 

91.5 


Se 

Ge 

2.341 

-44.10 

0.0344 

28.0 

Se 

Se 

2.330 

-40.60 

0.0318 

72.0 
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Se chain structure is maintained but the four fold tetra- 
hedially coordinated Ge atoms act as chain crossing points. 

In this model, Ge atoms are not allowed to bond to one 
another. The present results, obtained from the experimental 
phase functions, do not seem to support fully this particular 
moo cl . However it is clear that it is mostly Se which is 
found around Ge as well as Se and the probability of Ge 
Donding with another Ge is low. In the other model called 
the "random covalent model ( RCM ) " the four fold coordinated 
S' ■ atoms are allowed to bond with other Ge atoms as readily 
r r. the two fold coordinated Se atoms. The markedly differ- 
ing probabilities around Ge and Se do not warrant the vali- 
dity of RCM too. 

1 3 . i Proposed Model 

In determining the bond lengths from experimental 
■i ::;o functions it was assumed that a particular kind of 
.•■tom (either Ge or Se) can act as the backscattering atom, 
■■’M.. means that if Ge is four fold coordinated, either it 
V mended to all the four Ge atoms or to all the four Se 
atoms. However the possibility of both kinds of atoms 
being present in the first shell can not be excluded. For 
ex nmole Ge, which is tetrahedral ly coordinated to other 
four atoms, one can have five different possible configu- 
rations. The probabilities P Ge and P ge obtained from phase 
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ohift analysis themselves will determine the probability 
O'" -*- nc ^ v iclual configuration, i.e. if Ge is surrounded by 
four Ge atoms, probability of occurrence of this particular 
configuration would be P^. Figure 5.6 shows the different 
configurations for four fold Ge and two fold Se atoms with 
their respective probabilities. 

The probability of ith configuration around the 
central Ge atom is given by 


where M, 


1 

4 

6 


P . = 

l 

for i = 
for i = 
for i = 


M ± P 


i-1 

Ge 


1,5 

2,4 

3 


P 


5-i 

Se 


- 5 , 

U, p i = (P G S + p s e )4 = no 
i=l 


Then the first shell EXAFS \ (k) can be expressed 

5 

XjOc) = (5.4) 

i= 1 


where X^(k) represents the first shell EXAFS of ith 
configuration. 
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5 . 4 . 1 Calculation of Y (k) 

The expression for the first shell EXAFS is given 

by 

2r l 

N _ .22 

X T (k) = “■P 1 (k) e A e“ 2 ^l k sin f2kr 1 +4 > i (k 3 


',’hore is the number of particular kind of atoms surround- 
ing the central atom at an average distance r^. <|o(k) is 
the phase shift calculated by taking the backscattering 
nhase shift of that particular kind of atom. The expression 
can bo visualised as sum of terms y with each term repre- 
senting contribution from individual bond. Thus the first 
shell EXbPS for i -t h configuration can be calculated by 
taking such a sum. For example if a particular configura- 
tion has number of Ge neighbours equal to and number 

o.f. Je neighbours equal to N ge such that N Qe + N ge = 

then 


■/-,(k) = ^N Ge F Ge (k) sin [2kr x + ^(k)] + N ge Tg e (k) 


sin |2kr 1 + <b ge (kjj j- 


2 r 1 j P 

~2<Tk^ 

e 


kr-, 


where F„ (k) and F„ (k) are the backscattering amplitudes 

for Ge and 3e respectively, 4^ Ge ( k ) and ^Se^ are tKe 
total phase shifts calculated using the same central atom 
phase shift but different backscattering atom phase shifts 
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corresponding to Ge and Se respectively. Here the Ge-Ge 
and Ge-Ge bondlengths (or Se-Se and Se-Ge) have been 
assumed to be same and equal to r^. In evaluating equation 
5.4 it was assumed that the temperature dependent Debye- 
Waller type factor remains same for all the configurations 
and hence could be pulled out of the summation. The mean 
free path term could be absorbed in an overall scaling 
factor. 

5.4.2 Determination of structural parameters from the_model 

The first shell EXAFS was extracted from the experi- 
mental data as described in section 5.2 with E q = E c * The 
method of determination of structural parameters involves 
fitting the theoretically evaluated first shell EXAFS 
(equation 5.4) with the experimental first shell EXAFS with 
n set of four parameters namely the change in the threshold 
E , the first shell radius r ±/ the Debye-Waller type 
factor 2 and an overall scaling factor N. A least squares 
minimization routine MINFUN (see Appendix I ) was used 
for this purpose. In figures 5.7 and 5.8 are compared 
the experimental first shell EXAFS (associated with Ge edge 
and Se edge in II respectively) with the one obtained 
from the proposed model. The results are summarised in 
table 5.3. The average first shell radius around Ge as 
well as Se in both I and II turns out to be 2.39 ± 0.01 1 
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Table 5.3 

Parameters obtained from proposed model 


Compound 

Central 

atom AE Q (eV) 

„ 

2 °2 
£> ( A ) 

N 

I 

Ge 

37.5 

2.39 

0.0055 

-0.120 

I 

Se 

-4.3 

2.40 

0.0054 

-0.0688 

II 

Ge 

15.3 

2.39 

0.0044 

-0.099 

II 

Se 

-52.6 

2.39 

0.0040 

-0.058 
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The values of -£E e are consistent with their respective E 

2 

values. The factor w remains constant for particular 
compound irrespective of the kind of central atom. The 
overall scaling factors obtained for different cases could 
not be interpreted in simple way due to the various factors 
which affect the EXAFS amplitudes. 

b . 5 Results : DRDF 

X-ray diffraction patterns for both I and II have 
been recorded. Broad halos and absence of any sharp dif- 
fraction peaks indicated the amorphous state of -the samples. 
Figure 5.9 shows the diffraction pattern for II. The 
experimental data for DRDF analysis consisted of the dif- 
fracted intensity from the sample (with substrate), I, and 
from the substrate, I Q . The data were corrected for substrate 
scattering, absorption, polarization, compton scattering 
and dispersion of form factors. The data was then normalized 
to electron units by high angle method (Haider et al., 1966). 
The interference function F (s ) is obtained by equation 

f(s> = 3 \i co h - <f2> > < f>2 l 

where I is the first ordered coherently scattered inten- 
sity obtained by removing the modified component from the 
normalised intensity, 

< f > = I x i f i (s) ' <f2> = ^ x i f i (s) 

i 1 



i n 




at 



fig. $.10 Intensity in «l#ctron unit* for IX 
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with and f^(s) representing the atom fraction and 
scattering factor of element i respectively and 
g = 4/c'sine/x . 

The interference function F(s) is related to the 
radial distribution function as follows : 

F ( s ) = ( G(r) sin sr dr 

■J 

O 

where G(r) equals 4 ;a r( f - ? ). V is the average 
density and r is the radial distance. Fourier inversion 
of the above equation yields 

G(r) = — \ F(s) sin rs ds 

The termination of data at s max < co introduces 
spurious details in the resultant G(r). Following the 
p rocodur? for termination correction (Kaplow et al., 1965), 
th i final G(r) function was determined. 

To illustrate, the normalised intensity and the G(r) 
r unction Cor sample II are shown in figures 5.10 and 5.11 
respectively. The near neighbour distance from G(r) func~ 
Lion indicate a value 2.38 ± 0.03 1. Similar analysis 
Cor cample I also yields the same near neighbour separation. 



CHAPTER VI 


SUMMARY 
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The present work includes the structural study of 
crystalline and amorphous germanium and amorphous Ge Se. 

Jx _L 

system by X-ray diffraction. X-ray fluorescence and EXAFS 
).• •ehnigues. More emphasis is placed on the relatively new 
L 'otnique - EXAPS. Although the revival of interest in 
X/ i ihS is due to the improvements in theory and experimenta- 
tion, Ln particular the synchroton related developments, 

■;-> have in our laboratory recorded the EXAFS data with a 
*ttt< vantional X-ray source. By operating the tube at higher 
Lr t. In is the requirement of high photon flux was met. However 
this would introduce systematic errors in the experimental 
(is La us the harmonics contribute to the measured intensity. 

1 L has been assumed, while analysing the EXAPS data, that 
;.iis harmonic contribution is relatively small and thus 
could be neglected. Adequate resolution (~»4eV) in the 
vi-’inity of absorption edge was obtained by using an exit 
u r ; mm width (which allows minimum increment of 0.01° 

in Because of the slight variation in the intensity 

u ;■ radiation from the source with time, it is expected that 
;-t, P incident and transmitted intensities (I Q , I) be mea- 
cur-d simultaneously. However in the present work we have 
m insured the transmitted intensity (in preset count mode, 

1 0" 1 counts) at all the experimental points (441 points for 
Go idge and 331 points for Se edge) and then repeated the 
scan without the absorber in the path of the beam. The 
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variation in the intensity with time was found to lie well 


,.;il;!iin the statistical precision. 


The analysis of EXAPS data involved number of steps 
;ith expensive use of computer. The programmes needed for 
LI.'"' analysis were developed (see Appendix I ). The experi- 
mental data of crystalline and amorphous germanium served 
to standardize the analysis procedure. In a bid to determine 
the near neighbour distance by direct comparison of experi- 
mental and theoretical phase shifts it was found that it is 
possible only when the phase functions are extracted at 
LC = E When near neighbour distances were determined 
with values of E q other than E c , no physical significance 
could be attached to the parameter &E q( the change m the 
Threshold. The distances determined also varied by as much 
, s 0.01 A°. On the other hand analysis with E q = E c yields 
U.r- ond lengths to an accuracy of 0.01 A°. Also the values 
o' AL 0 were of the order of E c with change m sign. This 
ul -, ;i5 that the threshold which was shifted by E c , to obtain 
Id ..or phase shift, is again shifted towards E q = o thereby 
Introducing non linearity in the phase shift. 


Calculation of theoretical phase shifts require the 
know led to of both the central atom phase shift and back- 
ccaU -ring atom phase shift. These have been tabulated 
by too and Lee (1979). The phase shifts calculated by 
then; using the He xman-Ski liman wave functions were not used 
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cis variation in the backscattering phase shift with atomic 
number posed a problem in interpolation of the tabulated 
data. The figure 6.1 shows the backscattering phase shifts 
at k = 3.7795 A° ^ calculated using Clementi-Roetti wave 
functions and Herman-Skillman wave functions (dashed curve) 
as a function of the atomic number. In our analysis we 
found the phase shifts calculated using Clementi-Roetti 
wave functions suitable. 

In the Fourier filtering technique which has been 
used to separate out the first coordination shell contri- 
bution to EXAFS, normally a 'sufficiently' large window 
encompassing the first peak is used. A systematic study 
has been carried out to illustrate the effect of the width 
of the filtering window on the phase and amplitude func- 
tions md subsequently on the local structure parameters 
determined using these f motions. The near neighbour 
distance determined from phase functions for different 
widths of the filtering window was found to remain constant 
(see Table 4.1 ). This indicates that the phase part of 
first shell EXAFS is hardly affected by the width of the 
filtering window. Contrary to this the amplitude function 
was found to be very sensitive to the window function. 

For - 'narrow' window the amplitude function simply dies 
away almost exponentially as the photoelectron wave vector 
increases. With increasing width of the filtering window 
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the trend of the amplitude function approaches to that of 
the backscattering amplitudes tabulated by Tea and Lee 
(1979). In our analysis we have used an optimum width 
(1.25 A°) which retains the details in the amplitude 
function (see section 4.6). 

Concepts of chemical transferability for both the 
phase and amplitude functions have been used while analysing 
the crystalline and amorphous germanium data. Transferability 
of phase functions yielded the near neighbour distance for 
a-Ge which was in agreement with the one obtained by compa- 
rison of theoretical and experimental phase shifts. The 
DRDF obtained from X-ray diffraction studies (see sec. 4. 10) 
also indicate almost the same near neighbour distance. 

' n ransf erability of amplitude functions was used to determine 

the coordination number of a-Ge. The logarithmic ratio of 

2 

the amplitude functions when plotted against k yields a 
non linear curve (see figure 4.12). The determined 
value of the coordination number was found to be very 
sensitive to the range of k-values used for fitting the 
straight line. Thus the coordination number could be deter- 
mine only with a limited accuracy (4.49 ± 0.8). The other 
method which has been used for determining the coordination 
number of a-Ge involves determination of an overall s c aling 
factor for c-Ge amplitude function and then using the same 

This method also gave a value of 


for a— Ge amplitude. 



121 


coordination number slightly greater than four with as 
much an accuracy as that in case of transferability of 
amplitude functions. The results obtained for a-Ge are 
.m agreement with the earlier RDF and EXAFS studies except 
foT - coordination number which appears to be greater 

then four, a result which is hard to credit. The reason 
j. o i- obtaining higher value of coordination number also 
suggests that the mean free path in case of a— Ge is large 
no compared with that in c-Ge. Earlier studies (Temkin 
ot al., 1973, 1972; Shvchik and Paul, 1972; Richter and 
hr 'Idling, 1951) have predicted coordination number of 
a-Gc to be less than four. Our results when corrected for 
moan free path and the accuracy with which the coordination 
m. J^r can be calculated with EXAFS is taken into account 
no , l<j not contradict the earlier results. 

In the present work, the study of Ge Se 1 system 

IX 1 -"iX 

1 ' ; ,,ocr 3 imited to compounds with low concentration of 

Tormanium. The earlier studies with X-ray diffraction 

nrc'osed two possible models namely the CCM and RCM. In 

rooont times many authors have supported the former model 

ro, Go te, svstem with xb0.33 (Ball and Chamberlain, 1978; 
x 1 -x 

d onons dock, 1976). In the present work we have studied 
t;h-> two composition x = 0.1 and x = 0.3. X-ray fluorescence 
measurements were used to ascertain the compositions of these 
compounds. The DRBF obtained from X-ray diffraction studies 
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of these samoles indicate the first neighbour separation. 

n 

to be equal to 2.38 A . This would correspond^ all the 
three bonds Ge-Ge, Ge-Se and Se-Se. 

In order to see whether EXAFS can give some informa- 
tion regarding the local structure, the EXAES associated 
with both Ge as well as Se absorption discontinuities has 
been studied. While EXAES typically measures first shell 
distances with an accuracy of 0.01 A°, the resolution of 
different bond lengths is much poorer of the order of 0.2A°. 
With much work and in favourable circumstances the limit 
can be pushed down to ^ 0.10 A° (Shulman et al., 1978). 

For the Ge Se. system thus the distances Ge-Ge and Ge-Se 

yz JL 

(or Se-Se and Se-Ge around Se) coud not be resolved. There- 
fore the first shell phase function could be on account of 
the two possible bondings. Out of the different procedures 
suggested for determination of first shell radius we choose 
t,'. to method of comparison of theoretical and experimental 
phase shifts. While comparing the phase shifts again one 
requires tho knowledge of backscattering atom. Depending 
on tii „* quality of fit obtained with each of the atoms as 
buckscattering atom and their concentration in the compound, 
v/o talk in terms of the probability of the particular atom 
b ling the nearest neighbour of the central atom. This 
approach was found to be quite satisfactory but does not 
support any of the earlier models. Therefore we make use 
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of these probabilities in the proposed model which describes 
the first shell EXAFS in - terms of all the possible configu- 
ration* but with different probabilities. The comparison 
of the experimental first shell EXAFS with the one calculated 
from the probability model yields the parameters which are 
quite consistent as discussed in section 5.4.2. Thus it 
has been possible to determine the distribution of atoms 
(in terms of probabilities) surrounding the absorbing atom 
using the EXAFS phase functions. 

The method developed in the present work for GeSe 

system may well be extended to study systems like As 2 Se 3 

or As Se system and the multicomponent systems such as 
x 1—x 

GeSeTe and so on. It is desirable to record the EXAFS 
data at liquid nitrogen temperature and this could help in 
the accurate determination of disorder parameter (Gregor 

and Lytle, 1979). 

Although the structure in FT beyond the first peak 
is significant, the present study was limited to only first 
coordination shell. However with the determined probabilities 
p and p t it is possible to simulate a structure which 
would explain the details in the Fourier transform beyond the 


first shell 
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APPENDIX I 
PROGRAMMING NOTES 

1 • Elementary Analysis 

A fairly simple programme was constructed (PROG— I) 
which reads the times T, T q for preset number of counts 
with absorber in the path of the beam and without the 
absorber. The subroutine VICTOR gives the Victoreen 
coefficients A and B. The subroutine INTER has been used 
to interpolate the data onto a grid of egually spaced k 
values (Kaplow et al., 1965). The input data is stored in 
file opened with unit 2 3. The output is written on X. CDR 
(unit 25) which gives the contribution to absorption 
exclusively due to K shell as a function of photoelectron 
wave vector k. D. CDR file (unit 24) contains the absorp- 
tion h x as a function of energy at all the experimental 
points. V. CDR gives the extrapolated curve A A + 3 ^ • 

7 # Fourier Transform 

A single subroutine FOUR (PROG— II) was used to perform 
the Fourier Transforms and Inverse Fourier Transforms 
required in the analysis of EXAFS. A NAG LIBRARY program 
C06FCF gives the Fourier coefficient a n ’s and b n 's. 
Combination of C06FCF and C06GCF (for taking complex conju- 
gate) results in inverse transforming the data. DELR is 
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the spacing in r at which the coefficients are calculated. 
It is given by 


DELR = Tv / (N*DEL ) 

where DEL is the spacing in k-values at which input data 
in tabulated at N points. 

3 . Minimization : MINE UN 

The program was written by W.E. Humphrey at the 
Lawrence Radiation Laboratory (W.E. Humphrey, "A general 
minimizing routine - MINFUN", Programmer' s notes P.6, 
h. 7 . 62, UCRL, Berkeley). The program was modified for 
some of its input requirements. Purpose of the program 
is to minimize a function of several variables. It can 
Do used for the following two distinct purposes : 

(a) Finding the minimum of a function of n parameters as 
a mathematical problem. 

( 'b) Finding a fit for a set of experimental data to a 
given function by minimizing the sum of chi-squares 
or negative logarithm of maximum likehood function. 
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PROG - I 

C ELEMENTARY DATA PROCESSING PROGRAM FOR EXAFS 

C PROGRAM READS N1 AS TOTAL POINTS, N2 NO. POINTS 

C FOR VICTOREEN FIT 

C N3 IS THE POINT CORRESPONDING TO MINIMUM VALUE OF I< 

C SWITCH IF POSITIVE, ENERGY, LAMDA AND MUE ARE 

C PRINTED AND PAUSE 

C MM IS NO. OF POINTS LESS 1 REQUIRED FOR FOURIER 

C FILTERING 

C TTHM LAST VALUE OF 20 

C DATA IS IN DECREASING ORDER OF 20 

C NBAD IS NO. OF BAD DATA POINTS AND NB IS ARRAY CONTAINI1 

C CORRESPONDING LOCATIONS. 

DIMENSION ALAMDA(500) , E ( 500) , AK( 1500) , AMUE( 1500), AA(144 
, AA1 (300,2 ),NB(10),TO(500),T(500) , TTH (500) , AMUEl (500), 
V ( 5 0 0 ) , AMUE2 (500 ) , AKl ( 500 ) 

OPEN (UNIT=23,DEVICE='DSK' ,FILE=’CGE.CDR' ) 

OPEN (UNIT=24 / DEVICE= 1 DSK* ,FILE=* D. CDR* ) 

OPEN (UNIT=2 5;DEVICE=' D3K' , FILE= ' X. CDR' ) 

OPEN ( UNIT=2 6, DEVICES DSK',FILE=' V. CDR* ) 

DATA IN, INI, IN2, IN3/2 3, 24, 2 5, 2 6/ 

READ( IN, 1) AA 

1 FORMAT (72A1/36A1/ 36A1 ) 

RE AD ( IN, *) N1 , N2 , N3, SWITCH, TTHM 
MM=1 02 3 

PRINT2 , N1 , N2 , N3, SWITCH, MM, TTHM 
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2 FORMAT (1 OX, 316, 5X, FlO. 4, 5X, 16, 5X,FlO. 4) 

READ (IN,*) (T(I) ,1=1, Nl) 

READ( IN, * ) (TO ( I ) , 1=1 , Nl ) 

READ ( IN, * ) EO , NBAD, (NB( I ) , 1=1 , NBAD) 

DO 3 1=1 # N1 

AMUE1 (I)=ALOG(T(I)/TO(I) ) 

TTH (I ) =TTHM- ( 1-1 ) *0. 01 

ALAMDA( I ) ( I ) =4 . 02 62 *SIND (TTH ( I )/2 « 0 ) 

3 E(I) =12398. 52/ ALAMDA( I ) 

C REMOW BAD DATA POINTS 

IP (NBAD. EQ. 0 ) G0T04 
DO 5 1=1, NBAD 
NB1=NB(I) 

5 AMUEl (NB1 ) =AMUEl (NB-1+1 ) 

4 WRITE (IN, *)N1 

WRITE ( INI, 90 ) (E (I ) , AMUEl ( I ) , 1=1, Nl ) 

90 FORMAT (10F10. 3) 

IP (SWITCH) 6, 7, 6 

7 PRINT1, AA 

PRINT8, (I, TTH ( I ) , E ( I ) , AMUEl (I) , ALAMDAd ) , 1=1 ,N1) 

8 FORMAT ( 3 3X, ' S . NO . ' , 8X, ' TTH ' , 1 IX, ' E ' , 15X, 1 MUET ' /// 
1 (30X,I6, 4F15.4) ) 

PAUSE 

6 CONTINUE 
READ(IN,*)NE1,NE2 



135 


SLl=( AMUEl (Nfil+1)-AMUE1 (NE1 ) )/(E (NEl+1 ) -E (NEl ) ) 
do 2 0 X =nEl +1 , NE2 

SL2 = ( AMUEl ( 1+1 ) -AMUEl (I ) )/(E( 1+1} -E(I)) 
PRINT*,I,E(I),SLl,SL2 
TYPE* , I, E (I ) , SL1 , SL2 
IP (SL2.LT.SL1) G0T021 
2 0 SJ1=SL2 

21 E01=(E(I-1 )+E(l) )/2.0 

S0=E01 

EXO 9 1=1 ,N2 

XI = AL AMDA ( I ) **3 

X2=X1*ALAMDA(I) 

AAl (I,l)=Xl 
9 AA1(I / 2)=X2 

CALL VICTOR(AAl,AMUBl,N2,2,C,D) 

PRINT*, C, D 

C EXTRAPOLATION AND SUBTRACTION 

DO 100 1=1, N1 
Xl=ALAMDA( I ) **3 
X2=X1*ALAMDA(I) 

100 V ( I ) =C*X1 +D*X2 

DO 10 I=N3,Nl 
J=I-N3+1 

AMUE2 (J)=AMUE1 (I)-V(I) 

AK1 ( J ) =SQ RT ( 0 . 2 6 3* (E(I)-EO) ) 
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10 CONTINUE 
N13=H1-N3+1 

DEL = ( 16. 0*3. 1415 92 65/857.0) 

Nl 1 1 = ( AK1 ( Nl 3 ) -AK1 (1 ) )/DEL+l 

CALL I NTE R ( AK1 , AMUE2 , Nl 3 , AK, AMUE , Nl 1 1 , DEL ,11) 
PAUSE 

WRITE (IN2, *)N111, DEL, II, AK(Nlll) 

WRITE ( IN2 , 1 1 ) ( AK( I ) , AMUE ( I ) , 1=1, Nl 11 ) 

WRITE (2 6, * )N1 

WRITE (26, 90) (E(l) , V(I), 1=1, Nl ) 

11 FORMAT(10P10,7 ) 

STOP 


END 
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PROG - II 

SUBROUTINE FOUR(X, Y, N, NL,NH, 11,12/ DELR) 

DIMENSION X(8192) ,Y(8192 )/WORK(8l92 ),XX(1000) ,Rl(1000), 
1WW(1000) 

DIMENSION WR( 500 ) , WR1 (100) ,WW1 (100) 

COMMON CONT 

OPEN (UNIT=50, DEVICE=' DSK' ,FILE=' R. CDR' ) 

OPEN (UNIT=51 , DEVICE=* DSK* ,FILE=’ W. CDR' ) 

X IS AN ARRAY CONTAINING REAL PART Y IS THE ARRAY 
CONTAINING IMPARINARY PART OF THE DATA TO BE ANALYSED 
IFAIL=0?PI=3. 141592 65 
CALL C06FCF (X, Y, N, WORK, IF AIL ) 

IF I1=NL=0 THE NEGATIVE R-PART IS REMOVED AND I.F.T. 
COMPLEX FUNCTION IS RETURNED 

IF 12=2 SMOOTH FILTERING WINDOW IS USED BEFORE I.F.T. 

REAL FUNCTION IS RETURNED 

IF (II .EQ. 3) GOTO 4 

IF (ML, EQ. 0) GOTO 5 

MMM=8 . 0/DEL R+l 

DO 6 1=1, MMM 

R1(I) = (I-1)^DELR 

XX ( X ) =SURT (x( I) * >2 +Y( I ) **2 ) 

CONTINUE 

IF ( 12 . EQ. 1 ) G0TO7 



8888 

7 

8 


ft 8 8 

786 

2222 

10 


WRITE ( 50, *)MMM 

WRITE (50, 8888) (R1 (i) ,XX(l) ,1=1, HMM) 

FORMAT( 10F10. 4 ) 

TYPE 8 

FORMAT (5X, ' TYPE 1 IF NL NH ARE KNOWN ELSE 0‘) 

ACCEPT*, TR 

IF (TR.EQ.1.0) G0T09 

CONTI =0. 0 

DO 786 1=1,500 

IF ( CONTI. NE. 0.0) GO TO 888 

IF(XX(I+1) .GT.XX(I) ) GO TO 786 

TYPE 11, (Rl(l),X(l),Y(l), I, XX (I) ) 

C0NT1=1. 0 
GO TO 786 

IF(XX(I+1) .LT.XX(I) ) GO TO 786 
CON T1 =0.0 
CONTINUE 
TYPE 10 

FORMAT(5X, ' GIVE NO OF COEFFICIENTS TO BE TYPED 
IF ROM N1 TO N2 1 ) 

ACCEPT*- , MM1 , MM2 
IF (MM2.EQ. 1) GOTOIOOOO 

TYPE 11, (R1 (I),X(I),Y(I),I,XX(I),I=MM1,MM2) 


10000 CONTINUE 



11 FORMAT(5X,3Fl6.7,I6.7) 

TYPE 5555 

5 55 5 FORMAT (’ TYPE 1 TO CONTINUE, 0 TO GO BACK') 
ACCEPT*, KK 

IF ( KK. EQ. 0 ) GO TO 2222 
TYPE 454 

4 54 FORMAT (' DO YOU WANT TO RETURN : Y/N' ) 

ACCEPT 1, CONT 
1 FORMAT (A4 ) 

IF (CONT.EQ. r Y* ) RETURN 
9 TYPE 12 

1 2 FO RMAT ( 5X, 1 GIVE NL NH' ) 

ACCEPT*, NL,NH 

1 • IF (ML. EO. 0 ) G0T013 

IF (NL. EU. 1 ) G0T088 

% 

X(l)=0.0;Y(l)=CM0 
DO 14 1=1 , NL-1 
X ( I-! 1 ) =0. 0; Y( I-l-l ) =0,0 
X(N-I+1 ) =0. 0; Y(N-I+1 ) =0. 0 
i 4 CONTINUE 

IF(NH.E0. 0)GOT04 
DO DO 16 I=NH“1 , N/2-1 

X(I+1)=0.0; Y(I+1)=0 
X(N-I+1 ) =0. 0; Y(N-I+1 ) =0, 0 
16 CONTINUE 
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IF (II. EQ. 2) GOTO 22 
A C^IA, C06GCP ( Y,N, IFaIL ) 

CALL. CQ6PCE (X, Y, N, WORK, IP AIL) 

C ALL CO 6GCP ( Y, N, IF AIL ) 

RETURN 

.1 '•) DO 2 0 I— 1 , N/2—1 

y (N— I+l ) =0. 0 ? Y(N-I+1 ) =0. 0 

20 CONTINUE 
GO TO 4 

2 2 rvWAX=SEA ( XX, MM, NL, NH ) 

DO 21 I=NL— 1 , NH-1 
11=1-1-1 

771*0. 6« ( 1-COS (2 . *PI* ( Il-NH) / (NL-NH) ) ) 
VJW( H 2 -NL ) =WI*XMAX 
NR i,I 1-2 -NL ) =Rl (1+1) 

X( I+l ) =X( I+l ) vWI; X(N— I+l ) =X (N—I-t-1 ) *WI 
V ( I+l ) =Y( I+l ) *WI ; Y(N— I+l ) =E (N-I+l ) *WI 

21 CONTINUE 
m,N-NH-NL+l 

ML A= 1 0 0 ; D = (WR(NLN) -WR (!))/( NLA-1 ) 
CALL INTER (WR, WW , NLN, WR1 , WW1 , NLA, D, II ) 

WRITE (51, *)NLA 

'WRITE (51, v ) ( WR1 ( I ) , WW1 ( I ) , 1=1 , NLA) 

GOTO 4 


end 



141 


FUNCTION S1A( A, N, Nl , N2 ) 
DIMENSION A(N) 

XMAX=A(N1 ) 

DO 10 I=N1+1,N2 
±F(XMAX.GE.A(l) ) GO TO 10 
»V0C=A(I) 

3 0 CONTIilUE 
oEA=XHAX 
AST URN 


END 
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APPENDIX II 

SEPARATION of phase and amplitude functions 
T he expression for the first shell EXAFS is 
N 1 If l , .2,2 

V. (k) = ■“ e“ i e"’ - ' K sin 1 2^-:- -f'(k)j 

kr“ 

--- A(k) sin -'f(k) (i) 

v/here A(k) and"-' (k) are the amplitude and phase functions 
respectively. 

The Fourier transform has the ability to decompose 
a s i ne wave into a unique amplitude and phase function. 

Wo can write (i) in the form 


. k) = 


2i 


A(k) e 


-f'(k) 


1 _ 

2i 


-i|0c) 


Cii) 


0 


If y ^ Or ) is Fourier transformed to F(r) and /--j_ (k) 
is sufficiently bandlimited in r (guaranteed by filtering), 
then the tsrm e"*" ' in (ii) corresponds to positive r 
values and e _i • ^ term negative ones. If, before trans- 
forming F ( r) back to ^(k), the negative r part of F (r) 
is replaced by zero, then resulting inverse transform 

p ro duces 


1 

2i 


A(k) 


i '■j' (k) 


z(k) 


/ 
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The phase and amplitude functions are readily obtained 
as 

A(k) = 2 j Z (k) | 

n|.j(k) = arg i Z(k)| + 



